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4.BSIMI 

A molecular  beam  apparatus  mas  used  to  investigate  the  scattering 
of  air  and  molecules  from  surfaces  of  polished  low  carbon  steel,  etched 
low  carbon  steel,  polished  aluminum,  unpolished  window  glass,  and  polished 
window  glass.  A beam  of  molecules  was  directed  against  the  surface  at  rep- 
resentative angles  of  incidence.  Polar  flux  distributions  of  the  scattered 
molecules  were  measured  by  a movable  detector.  An  ionization  gage  beam  de- 
tection system  was  devised  to  accomplish  this  measurement. 

Spatial  polar  plots  of  these  measured  flux  distributions  were  found 
to  approximate  closely  the  form  of  the  cosine  scattering  distribution  for 
the  cases  of  the  steel  and  aluminum  surfaces.  Similar  plots  in  the  case  of 
the  glass  surfaces  were  found  to  have  well  defined  but  small  bulging  devi- 
ations from  the  cosine  scattering  shape  on  the  side  away  from  the  incident 
beam. 

In  order  to  illustrate  one  aerodynamic  application  of  scattering 
data  obtained  by  molecular  beam  techniques,  a procedure  was  devised  for  the 
computation  of  the  coefficient  of  specular  reflection,  f,  (f  being  that  frac- 
tion of  the  total  incident  tangential  momentum  which  is  transferred  to  a 
surface  in  interaction  with  a flowing  gas)  and  values  of  this  coefficient  were 
calculated  for  the  case  of  an  unpolished  glass  surface  oriented  tangentially 
to  the  direction  of  mass  motion  in  a free  molecule  flow.  These  values  were 
calculated  as  a function  of  the  parameter  s,  the  molecular  speed  ratio  ( the 
ratio  of  the  gas  mass  speed  to  the  most  probable  thermal  speed  of  the  gas  mole- 
cules). If  all  scattered  molecules  are  assumed  to  have  the  distribution  in 
speed  in  a particular  direction  i&ich  they  would  have  had  had  they  issued  from 
a gas  in  equilibrium  at  the  temperature  of  the  surface,  f takes  on  values  of 
^ 0.97  as  s approaches  0 to  0.99  at  8 3 10.  Alternatively,  if  those 
molecules  contributing  to  the  deviations  from  cosine  scattering  are  assumed  to 
issue  from  the  surface  with  their  incident  speed  while  the  remainder  are 
assumed  to  be  diffusely  scattered,  (a  modified  Maxwell  model),  f takes  on  values 
of  'v  0.97  as  s approaches  0 to  0.93  at  s * 10,  These  values  do  not 
correspond  closely  to  the  published  values  for  air  on  glass  of  0.89  (svO), 
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tube  or  orifice  area.  (Subscripts  1,  2,  3,  etc.  de- 
note association  of  Cl  with  regions  of  the  beam  as 
required  in  section  2.31  and  2.32.  Subscripts  a,  b, 
etc.,  denote  the  a or  b end  of  a certain  tube  as 
required  in  sections  2.31  and  2.32). 
aS 

constant  ■ (Subscripts  r and  f,  if  used,  relate 

A to  reflected  molecules  or  to  molecules  issuing 
from  a fictitious  gas  respectively  as  described  in 
section  4.2). 

reflected  beam  signal  amplitude  in  arbitraiy  unitB 
for  incident  molecules  having  direction  ©{and  de- 
tector axis  having  direction  © # 

constant  representing  the  diameter  of  the  assumed 
cosine  scattering  shell. 

beam  signal  amplitude  differences  defined  ly 

-ADco$<t>' 

radii  of  tubes  or  orifices 
s sin  ©( 

magnitude  of  the  total  molecular  velocity  (molecular 
speed)  (subscripts  i and  r denote  incident  and  re- 
flected molecules  respectively.  Subscript  f denotes 
molecules  issuing  from  a fictitious  gas  as  described 
in  section  4.2). 


components  of  the  total  molecular  velocity. 


mean  molecular  speed  referred  to  coordinate  system 

fixed  to  the  flowing  gas  ~ Q, 

(Subscripts  , 1,  2,  3,  etc.,  denote  specific  mean 
speeds  as  required  in  section  2.31  and  2.32). 


mean  molecular  speed  of  molecules  issuing  diffusely 
from  a surface  - Z. 

4 ftr 

molecular  energy  fluxes  incident,  reflected  and 
diffusely  reflected. 


Maxwellian  distribution  function,  Eq.  4.3. 

coefficient  of  specular  reflection,  section  1.2  and 
Eq.  4*1 


magnitude  of  the  net  tangential  momentim  incident  on 
unit  surface  per  unit  time. 
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magnitude  of  the  net  tangential  momentum  incident  on 
unit  surface  per  unit  time  and  unit  solid  angle 
having  direction  , Eq.  4.25. 

magnitude  of  the  net  tangential  momentum  imparted  by 
the  surface  to  the  Issuing  molecules  per  unit  area 
per  unit  time,  Eq.  4.1  (the  addition  of  the  subscript 
2 indicates  the  application  of  assumption  E-2  in  the 
calculation  of  Gr). 

issuing  momentum  function  defined  Eq.  4*15. 

function  defined  in  Eq.  4*10. 

modified  Bessel  function  of  zero  and  first  order. 

Boltzmann's  constant 

molecular  mass 

number  of  beam  molecules  entering  gage  per .unit  time. 

number  of  molecules  of  all  velocities  striking  unit 
surface  area  per  unit  time,  Eq.  4.2. 

number  of  incident  molecules  per  unit  speed  range  and 
solid  angle  having  the  speed  and  incident  di- 
rection et,<h  striking  unit  surface  area  per 
unit  time. 

number  of  molecules  per  unit  solid  angle  and  unit 
speed  range  crossing  a unit  plane  in  a gas  in  a parti- 
cular direction. 

molecule  number  density  (the  absence  of  a subscript 
denotes  the  free-stream  density.  Subscripts  1,  2,  3, 
etc.,  denote  specified  regions.  Subscript  f denotes 
number  density  of  a fictitious  gas). 

instantaneous  value  of  molecule  number  density  in 
gage  at  time  t,  Eq.  2.321. 

distribution  function  denoting  the  probability  per 
unit  solid  angle  and  speed  range  that  a molecule 
with  incident  speed  ^ and  direction  0(|(pv  will  be 

scattered  with  speed  Cr  in  direction  ©r>r  (If 

present,  the  subscript  D identifies  the  probability 
distribution  function  with  the  diffusely  scattered 
molecules  while  the  subscript  S identifies  that  as- 
sociated with  the  specularly  scattered  molecules). 
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distribution  function  denoting  the  probability  per 
unit  solid  angle  that  a molecule  with  incident  di- 
rection ©i  !<fei,  unspecified)  will  be  scattered  in 
direction  0*.  (the  subscripts  D and  S perform 
the  same  function  as  above). 


distribution  function  denoting  the  probability  per 

unit  speed  range  that  a molecule  will  issue  from  a 

surface  with  the  speed  C -• 

r 

pressure  and  pressure  increment  (subscripts  if  used 
denote  specified  regions). 

surface  - detector  separation  (Fig.  7) 

beam  path  length  (Fig.  7) 

function  defined  Eq.  4.27 


molecular  speed  ratio 


0 

~Vwv 


free-stream  temperature,  Kelvin  (subscripts  1,  2,  3, 
etc.  if  applied  refer  to  specified  regions.  Sub- 
script f refers  to  fictitious  gas.  Subscript  w re- 
fers to  wall  temperature). 


magnitude  of  the  mass  velocity- 
detector  gage  volume 
molecular  most  probable  speed  = 


]T^1  _ 

1Z% T*  - 


I 

15 


(Subscripts  r,  f,  if  used,  associate  7m  with  T^, 
respectively) . 


probability  that  a molecule  striking  the  plane  of  the 
end  of  a tube  of  specified  shape  from  a gas  at  rest 
will  pass  through  without  returning  to  the  entrance 
plane  (subscripts  1,2,3,  etc.,  when  used,  denote 
specific  tubes). 


coordinate  system  of  incident  momentum  considerations. 

Fig.  41 

coordinate  system  of  issuing  momentum  considerations. 
Fig.  42. 
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accommodation  coefficient,  section  1.2 
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(subscripts  perform  the  same  function  as  for  V,wv ) 


incident  and  reflected  molecule  direction  in  unprimed 
coordinate  system , Fig.  41* 


reflected  molecule  directions  in  primed  coordinate 
system , section  4 .4. 


X 


mean  free  path 


A 

V 


microns  Hg  = 0.001  mm  Hg,  a measure  of  pressure, 
coordinate  system  fixed  in  a gas. 

polar  angle  as  measured  from  axis  in  above  coordi- 
nate system. 

element  of  solid  angle  in  above  coordinate  system. 


dto,  <Ato‘ 


element  of  solid  angle  in  unprimed  and  primed  X,Y,Z, 
coordinate  systems. 


SYMBOLS  USED  IN  SECTIONS  2.33  and  2.34  ONLY 


A = 
Ei  * 

Eo  = 
e,  * 

6 * 

"R.  3 2.  etc.  = 
T,  etc.  = 
V, , Vt  etc,  = 


amplifier  loop  gain 

potential  across  control  resistor,  emission  regu- 
lator. 

reference  potential,  emission  regulator, 
error  signal,  emission  regulator 
signal  to  meter,  difference  amplifier, 
signal  current,  difference  amplifier, 
circuit  resistances, 
circuit  transformers, 
circuit  vacuum  tubes. 
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AN  EXPERIMENTAL  MOLECULAR  BEAM  INVESTIGATION 
OF  THE  SCATTERING  OF  MOLECULES  FRCM  SURFACES 

1.0  INTRODUCTION 

1.1  General 

Since  the  time  of  Maxwell  and  extending  through  the  first 
third  of  this  century  there  have  been  many  investigation  in  the  field 
of  rarefied  gases.  In  recent  times  there  has  been  a resurgence  of 
interest,  stimulated  in  part  by  the  need  for  a wider  tinderstanding 
of  aerodynamic  conditions  met  by  bodies  in  high  speed  travel  through 
the  upper  atmosphere.  Of  the  many  investigations  in  this  field  only 
a small  portion  have  been  concerned  directly  with  the  interaction  of 
gas  and  surface.  It  was  the  objective  of  the  program  of  this  report 
to  investigate  one  aspect  of  these  gas-3urface  interactions  by  mole- 
cular beam  technique. 

1.2  Low  Density  Flow  Investigation  as  a Technique  of  Gas-Surface 

Inters  rtten  Stedy 

Historically  the  problems  of  the  interaction  of  gas  and  sur- 
face were  first  investigated  following  the  realization  by  Kundt  and 
Warburg  that  the  gas  was  not  at  rest  at  the  boundaries  of  rarefied 
gas  flows  (Ref.  l).  This  mass  motion  is  determined,  in  part,  by  the 
detailed  behavior  of  the  molecules  in  interaction  with  the  surface. 
Considerations  relating  to  these  phenomena,  therefore,  have  been  of 
great  interest  to  aerodynamicists  and  have  played  an  important  role 
in  the  formulation  of  rarefied  gas  flow  boundary  conditions.  Con- 
versely, from  experimental  investigations  of  rarefied  gas  flows  for  which 
model  boundary  conditions  have  been  formulated,  some  deductive  under- 
standing of  the  gas  surface  interactions  has  been  gained.  Summaries 
of  the  major  concepts  and  results  of  rarefied  gas  flow  theory  and  ex- 
periment are  contained  in  references  2 and  3« 

Maxwell  in  1879,  (Ref.  4)  proposed  a hypothetical  model  for  the 
gas-surface  interaction  with  respect  to  momentum  transfer  which  was 
consistent  with  the  available  data  and  which  has  been  used  in  part  in 
the  formulation  of  boundary  conditions  to  this  date.  If  some  molecules 
are  assumed  to  be  scattered  diffusely  from  the  surface  at  some  tempera- 
ture as  yet  unspecified  such  that  the  scattered  flux  distribution  is 
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in  accordance  with  Lamberts  law  of  diffuse  scattering  without  re- 
gard to  incident  angle  or  speed,  and  if  the  remainder  are  supposed 
to  be  reflected  specularly,  then  one  can  define  f , termed  the 
■coefficient  of  specular  reflection*,  as  that  fraction  of  all  mole- 
cules which  is  diffusely  scattered.  Since  only  this  fraction  contri- 
butes to  the  tangential  momentum  transfer  to  the  surface,  one  is  led 
immediately  to  an  alternative  definition  of  f , 

f - Gl  - Gr 

Gt 

where  Gi.  is  the  total  incident  tangential  momentum  broughtto  the 
surface  by  all  molecules  and  Gy  is  that  carried  from  the  surface. 

Only  those  molecules  which  are  specularly  reflected  in  the  Maxwell 
model  contribute  to  Or  . The  coefficient  f can  be  restated,  then, 
as  being  that  fraction  of  the  total  incident  tangential  momentum  which 
is  transferred  to  the  surface.  It  is  in  this  sense  that  f is  most 
reasonably  defined  since  it  can  be  applied  to  a macroscopic  experiment 
without  the  making  of  assumptions  regarding  the  detailed  interaction. 
It  is  this  definition  which  will  be  employed  in  this  report. 

In  its  past  application  f has  been  assumed  to  be  constant  of  the 
surface  and  the  gas  independent  of  gas  mass  velocity  or  gas  and  sur- 
face temperature  in  the  absence  of  more  comprehensive  data. 

As  a result  of  somewhat  analogous  reasoning  applied  to  the  de- 
velopment of  a theory  of  heat  conduction  in  gases,  Knudsen  in  1911 
(Ref.  5)  proposed  a coefficient  o<  , the  accommodation  coefficient, 
defined  as 

« m Ei  - Er 

El  - Ew 

where  E{,  is  that  energy  brought  up  to  the  surface  by  the  incident 
molecules,  E^is  that  which  is  carried  away  by  them  as  they  are 
scattered  from  the  surface,  and  is  that  energy  which  would  be 
carried  away  by  the  issuing  molecules  if  they  had  come  from  a gas  in 
Maxwellian  equilibrium  at  the  temperature  of  the  wall.  The  accommo- 
dation coefficient  has  been  in  use  to  the  present  in  the  formulation 
of  boundary  conditions.  As  with  £ it  has  been  assumed  to  be  constant 
of  the  surface  and  the  gas. 
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There  exists  a small  body  of  literature  pertaining  to  theo- 
retical and  experimental  determinations  of  f as  a result  of  studies 
of  macroscopic  flows.  A table  of  all  known  values  of  f from  such 
determinations  may  be  found  in  Table  I. 

Although  the  work  of  this  report  is  applicable  only  in  considera- 
tion of  momentun  transfer,  a table  of  experimental  values  of  ot  is 
included  in  Table  II.  Many  important  investigators  whose  names  do 
not  appear  in  the  body  of  this  report  are  referenced  in  these  tables 
and  so  appear  in  the  bibliography. 

It  can  readily  be  seen  from  these  tables  that  our  knowledge  of 
either  of  these  quantitites  is  not  extensive.  More  objectionable 
than  this  is  the  likelihood  that  the  values  as  reported  are  valid 
only  for  flow  conditions  approximating  those  of  the  experimental 
determinations.  The  chief  deficiency  of  these  determinations,  how- 
ever, lies  in  the  fact  that  they  yield  veiy  little  Insight  into  the 
physical  nature  of  the  gas-surface  Interaction  beyond  the  essential 
fact  that  by  and  large  the  molecules  "forget"  their  past  history 
upon  interaction  with  a surface. 

1.3  Detailed  Investigations  of  Gas -Surface  Interactions 

Soon  after  the  first  concepts  of  the  gas-surface  interactions 
were  formulated  investigations  were  undertaken  to  study  these 
phenomena  in  more  detail.  Some  variation  of  the  molecular  beam 
technique  was  almost  invariably  used.  These  investigations,  that  is, 
involved  a directed  stream  of  molecules  in  interaction  with  a plane 
surface  element.  Measurements  of  the  flux  of  scattered  molecules 
at  one  or  more  positions  about  the  scattering  surface  were  obtained. 

A considerable  body  of  literature  exists  devoted  to  the  description 
of  these  experiments  and  to  theories  of  the  interaction.  Two  very 
comprehensive  surveys  of  the  techniques  and  results  of  experiment 
are  contained  in  Befs.  2,  and  60  A summary  table  of  the  best  known 
of  these  investigations  can  be  found  in  Table  III.  Those  investi- 
gators not  specifically  mentioned  in  the  body  of  the  report  are 
referenced  in  that  table.  Beyond  pointing  out  certain  salient  fea- 
tures of  these  investigations  it  would  not  further  this  report  to 
elaborate  on  them  in  greater  detail. 
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Of  these  many  investigations  the  most  definitive  was  clearly 
that  of  I.  Estermann  and  0.  Stem  (Ref.  7)  in  their  classical  ex* 
perlments  on  the  diffraction  of  He  and  H2  from  cleavage  planes  of 
LiFl  and  NaCl.  The  location  of  the  observed  diffraction  maxima  mas 
found  to  be  in  strict  accordance  with  the  prediction  of  the  optics 
of  cross  gratings. 

The  majority  of  the  remainder  of  the  investigators  also  worked 
with  cleaved  surfaces  of  natural  crystals.  For  the  most  part  the 
mode  of  molecular  beam  detection  dictated  the  use  of  various  metal- 
lic atoms,  more  or  less  readily  condensible  at  the  temperature  of 
the  scattering  surface,  for  the  parent  beam  material.  These  experi- 
ments did  not  in  general  permit  such  a clear  and  unambiguous  inter- 
pretation of  results  as  did  the  work  of  Estermann  and  Stem. 

Only  one  group  of  tests  on  surfaces  other  than  crystal  cleavage 
planes  appears  appropriate  to  this  present  report.  These  were  per- 
formed by  Knauer  and  Stem  (Ref.  8)  on  polished  speculum  metal. 

An  optical  criterion  for  the  specular  reflection  of  transverse 
waves  is  that  the  irregularities  of  the  reflecting  surface  be  less 
than  a wave-length  in  extent  as  measured  normal  to  the  axis  of  wave 
propagation.  This  condition  was  believed  to  occur  at  glancing  angles 
of  0.001  radian  and  less  for  Hg  and  the  polished  metal.  Experiment 
verified  these  predictions.  Above  0.001  radian  no  detectable  amount 
of  specular  reflection  was  observed. 

On  the  basis  of  the  foregoing  remarks  and  Table  III  the  state 
of  knowledge  of  the  surface  interaction  may  be  summarized  as  follows! 

(1)  He  and  Hg  interact  with  cleavage  planes  of  simple  crystals 
in  a fashion  predicted  by  wave  mechanics. 

(2)  Gases  having  high  forces  of  interaction  with  the  crystal  sur- 
face apparently  scatter  diffusely  from  it.  For  example,  Li, 

Na,  Cs,  from  NaCl. 

(3)  Apparently  there  is  a tendency  for  some  metallic  gases  to  undergo 
a classical  elastic  specular  reflection  at  a crystal  surface 
where  the  interaction  forces  are  low  in  a manner  not  readily 
understood  in  wave  mechanical  terms. 

(4)  Gas-surface  Interactions  involving  surfaces  which  are  “liquid" , 
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micro -crystalline  or  man-formed  by  polishing  have  been  studied 
hardly  at  all*  The  past  findings  in  this  area  suggest  that  the 
scattering  is  largely  diffuse  at  least  wherever  the  criterion 
for  optical  specular  reflection  is  not  satisfied* 

There  la,  however,  experimental  evidence  from  the  determination 
of  values  of  f that  some  interactions  are  not  entirely  diffuse* 

For  example,  Millikan  reported  a value  of  f ■ 0,89  for  air  on  glass* 
Clearly  it  is  here  in  these  areas  of  diffuse  or  semif -diffuse  scat- 
tering that  our  understanding  of  the  molecule-surface  interaction 
is  most  meager* 

1*4  Objectives  of  the  Experiment 

In  most  general  terms  it  is  the  objective  of  the  molecular  beam 
program  to  study  the  fundamental  nature  of  the  gas-surface  inter- 
actions by  molecular  beam  techniques*  More  specifically  the  present 
experiment  was  designed  to  enable  the  measurement  of  the  issuing  flux 
distributions  of  atmospheric  constituent  molecules  scattered  from  sur- 
faces of  certain  ordinary  materials  oriented  at  arbitrary  angles  to 
the  incident  molecular  beam* 

It  was  recognized  that  the  work  herein  reported  could  constitute 
only  a beginning  to  a systematic  inquiry  into  the  nature  of  the  mole- 
cule-surface interaction*  Therefore* in  establishing  the  choice  of 
gas  and  surface,  experiment  geometry  and  resolution,  an  attempt  was 
made  to  be  guided  by  consideration  of  what  should  constitute  a logi- 
cal first  step  in  such  an  investigation.  The  details  of  the  experi- 
mental technique  and  of  the  test  program  will  be  dealt  with  in  the 
appropriate  sections  of  this  report. 

In  addition  to  reporting  the  results  of  the  experimental  pro- 
gram and  discussing  these  findings,  it  was  a further  objective 
of  this  paper  to  suggest  by  analysis  the  desirability  of  restating 
the  rarefied  gas  flow  boundary  condition  for  momentum  transfer  in 
terms  of  a more  general  transfer  ratio  than  the  constant  f * In 
connection  with  this  analysis  a technique  was  developed  for  computing 
the  value  of  f , subject  to  certain  assumptions  to  be  discussed 
later  in  this  report,  for  a Maxwellian  gas  with  superimposed  mass 
velocity  tangentially  directed  relative  to  a surface.  This  compu- 
tation was  carried  through  for  the  N^-glass  interaction  on  the  basis 
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of  data  obtained  aa  a portion  of  the  experimental  work  of  this  re- 
port. 

2.0  INSCRIPTION  OF  THE  APPARATUS 

2.1  Iq-fcreflugtaiy  Baaifcg 

The  majority  of  molecular  beam  devices  in  existence,  and  those 
most  familiar  to  physicists  today,  have  been  constructed  for  the 
study  of  the  fundamental  properties  of  the  atom  or  nucleus  as  re- 
vealed by  the  interaction  of  these  particles  with  electromagnetic 
fields.  Since  it  might  be  expected  that  a molecular  beam  apparatus 
constructed  for  the  study  of  the  scattering  of  molecules  from  sur- 
faces would  differ  slightly  from  these  other  equipments,  a brief 
description  of  the  equipment  of  this  program  will  be  given. 

The  molectilar  beam  is  formed  by  the  thermal  effusion  of  mole- 
cules through  an  orifice  or  tube  from  a small  source  chamber.  The 
effusing  molecules  travel  rectilinear  diverging  paths  to  the  plane 
of  the  defining  orifice.  Those  traveling  in  the  proper  direction 
pass  on  through  and  constitute  the  molecular  beam.  The  remainder 
(the  vast  majority)  are  stopped  and  become  removed  by  the  source 
pump. 

The  molecular  beam  downstream  from  the  defining  orifice  travels 
through  a region  of  high  vacuum  (P'V  10”^  mm  Hgj  mean  free  path, 

50  meters)  until  it  strikes  the  scattering  surface  oriented  at  an 
externally  variable  angle  of  incidence  to  the  beam  axis.  Those 
molecules  issuing  from  the  surface  into  a small  solid  angle  A CO* 
at  O'  and  (J)1  measured  from  coordinate  axes  fixed  in  the  surface 
(Figure  l)  pass  into  the  detector  and  create  a beam  signal  by  Pro- 
ducing a small  pressure  increment  within  that  gage.  By  moving  the 
detector  through  all  values  of  the  coordinates  one  can  obtain  the 
relative  distribution  of  molecular  flux  into  all  elements  of  *the 
hemisphere  at  a fixed  value  of  the  incident  angle.  By  repeating  the 
measurements  for  representative  values  of  the  incident  angle  a des- 
cription of  the  scattering  of  the  molecules  of  a given  gas  from  a 
given  surface  can  be  obtained.  The  data  of  this  report  were  obtained 
in  this  fashion.  Note,  however,  that  this  description  gives  a pic- 
ture of  scattering  behavior  which  is  by  no  means  complete,  since  it 
applies  to  one  incident  gas  temperature  only,  one  surface  temperature 
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only,  and  makes  no  measurement  of  the  distribution  in  speed  of  the 
scattered  molecules* 

A description  of  the  molecular  beam  apparatus  is  most  logically 
divided  into  two  principle  sections,  one  concerning  the  physical 
apparatus,  and  the  other  concerning  the  molecular  beam  detection 
system, 

2.2  Tfte.  PtoBlral  ipparatag 

In  briefest  terms  the  physical  apparatus  consisted  of  the 
beam  forming  and  detecting  elements  supported  within  vaeum  en- 
velopes which  were  in  turn  associated  with  an  evacuating  system* 

Figure  2 is  a schematic  of  the  pumping  system  showing  the  source  and 
detector  regions  of  the  apparatus  and  attached  pumps,  valves  and  re- 
frigeration components*  Figure  3 is  a photograph  of  the  complete 
equipment.  The  principle  vacuum  envelopes  (source  region  and  de- 
tector region)  were  of  stainless  steel.  Other  parts  were  variously 
of  cold  rolled  steel,  aluminum,  copper,  or  brass.  Nearly  all  gas- 
keting was  by  means  of  rubber  *0“  rings  for  fixed  as  well  as  ro- 
tating seals.  With  the  exception  of  the  quarter  swing  valves  be- 
tween the  main  envelopes  and  the  refrigerated  traps,  all  valves  were 
equipped  with  Teflon  seats*  Under  beam-on  conditions  typical  de- 
tector chamber  pressures  were  of  the  order  of  10*"^  mm  Hg  or  less 

-5 

while  the  source  chamber  pressures  were  of  the  order  of  2 x 10  mm 
Hg  or  greater. 

All  parts  of  the  molecular  beam  forming  and  detecting  equipment 
were  designed  to  be  supported  on  an  optical  bench  which  in  turn  was 
supported  by  the  plate  dividing  the  source  vacuum  envelope  from  the 
detector  envelope.  The  mounted  elements  were,  in  order  of  their 
appearance  in  the  beam  forming  or  detecting  process,  the  source 
chamber  and  source  orifice,  the  defining  orifice  (carried  by  the  bench 
support  plate),  the  shutter,  the  secondary  defining  orifice,  the  de- 
tector mechanism  (with  test  surface  and  detector  gage),  and  the  com- 
parator gage.  Attached  to  the  downstream  end  of  the  optical  bench 
was  a plate  (the  "apron")  carrying  coordinate  counters  and  screw 
driver  tipped  fittings  linking  the  moveable  parts  on  the  optical  bench 
with  the  rotating  shafts  piercing  the  detector  envelope  and  plate* 
Figure  4 shows  the  optical  bench  and  all  mounted  components* 
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The  source  chamber  mas  mounted  on  a slide  providing  adjustment 
of  the  source  orifice  center  by  means  of  lead  screws  and  ways.  The 
source  "orifice"  was  a tube  mounted  so  as  to  make  possible  precise 
allignment  of  the  beam  axis  with  the  source  tube  axis.  A pair  of 
auxiliary  optical  bench  carriages  with  transverse  and  vertical  slides 
and  coordinate  counters  and  mounting  a Point-o-lite  and  a telemicro- 
scope respectively  were  used  for  the  beam  alignment  throughout. 

The  source  chamber  was  connected  to  the  source  gas  supply  system  and 
manometer  system  shown  in  Figure  5. 

The  most  complex  of  the  various  components  of  the  beam  forming 
and  detecting  system  was  the  detector  mechanism.  Basically,  this 
mechanism  was  a mounting  device  for  holding  two  elements  of  the  mole- 
cular beam  system!  i.e.,  test  surface,  and  detector  gage,  in  the  de- 
sired mutual  geometrical  relationship.  The  test  surface  was  mounted 
normal  to  the  axis  of  the  aximuth  rotation  axis  by  means  of  a shaft 
along  that  axis.  (See  Fig.  6).  The  test  surface  and  all  other  ele- 
ments above  the  transverse  carriage  could  be  rotated  in  a horizontal 
plane  about  a vertical  line  element  in  the  test  surface  face  by  means 
of  exterior  controls.  Thus  the  incident  angle  of  the  beam  with  res- 
pect to  the  test  surface  was  variable.  The  upper  elements  of  the 
detector  mechanism  supported  the  detector  gage.  These  elements  were 
rotatably  mounted  to  provide  adjustment  of  the  polar  and  azimuth 
angles  of  the  gage  entrance  tube  referred  to  the  polar  axis  normal 
to  the  test  surface.  Azimuth  angle  control  was  made  by  internal  ad- 
justment. The  polar  angle  was  controlled  externally  by  means  of 
flexible  shafting  and  swiveling  miter  gear  linkages.  The  gage  sup- 
port system  was  terminated  by  a small  ground  bar  way  which  carried 
the  gage  holder.  The  gage  holder  provided  adjustment  of  gage-to- 
surfaee  separation  and  rotational  and  vertical  allignment.  All 
motions  not  considered  "set-up"  were  provided  with  coordinate  indi- 
cators in  the  form  of  5 place  Veeder  Root  counters.  Angular  position 
settings  were  determined  to  have  a reproducibility  of  0.015* o On  the 
basis  of  guaranteed  minimum  gear  quality  the  total  cumulative  error 
in  angular  setting  was  estimated  to  be  approximately.  0.05*.  For  the 
incident  and  polar  angle  gear  trains  the  counters  rotated  38.519  units 
per  degree. 
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Directly  associated  with  the  physical  apparatus  sere  certain 
Items  of  general  instrumentation  Including  pressure  monitors,  source 
gas  regulation  and  pressure  measurement  devices,  and  protective  de- 
vices. Three  Piranl  type  vacuum  gages  were  employed  to  monitor  pres- 
sures in  the  appropriate  range,  one  being  connected  to  the  fore 
pressure  line,  one  between  the  booster  pump  and  the  diffusion  pumps, 
and  one  to  the  source  pressure  manifold.  Both  the  source  and  detector 
sections  of  the  main  vacuun  envelope  were  provided  with  Alpert  type 
ionization  gages  for  base  pressure  monitoring.  Each  was  connected 
to  a standard  Low  Pressures  Project  design  power  supply  and  metering 
circuit.  This  device  is  provided  with  an  over-pressure  protection 
circuit  which  removes  filament  power  if  the  pressure  should 
exceed  the  maximum  pressure  of  the  range  selected.  At  the  same 
time  power  is  removed  from  a control  line  which  may  be  connected  to 
other  devices.  In  this  way  the  detector  and  comparator  gages  were 
also  protected  against  the  effects  of  over-pressure  conditions. 

The  source  pressure  manometer  indicated  in  Figure  5 was  an 
early  forerunner  of  the  present  Low  Pressures  design  precision  oil 
manometer  (Ref.  25).  While  it  lacked  certain  mechanical  and  optical 
refinements  of  the  present  models,  oil  column  height  differences  were 
easily  read  to  within  0.001  inches.  The  reliability  of  readings  in 
calibration  experiments  with  a precision  McLeod  have  been  found  to 
be  - 1 micron.  A more  detailed  discussion  of  the  design  consideration 
leading  to  the  present  equipment  may  be  found  in  Ref.  36. 

2.3  IhQ.  Jol9S]ter.j9flfl  J>9.te9,U<2B  J&fl&gtt 
2 «3i 

A portion  of  the  molecular  beam  apparatus  necessitated  by 
the  present  experimental  program  was  a beam  detection  system 
for  atmospheric  constituent  molecules.  Such  a system  must  in- 
clude a gage  which  collects  a number  of  directed  molecules  and 
measures  the  resultant  increment  in  molecule  number  density 
within  its  volume. 

Referring  to  Figure  7,  with  the  gage  set  to  receive  the 
direct  beam  and  with  the  test  surface  withdrawn,  let  the  mole- 
cular beam  shutter  be  opened  suddenly,  the  molecule  number 
density  within  the  gage  having  been  at  equilibrium  with  the 
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residual  gas  within  the  vacum  shell • ThenN  directed  mole- 
cules per  second  pass  into  the  gage  volume.  At  equilibrium 
with  this  new  influx  ^3  C-3  03b  molecules  per  second  strike 
the  area  03b  from  within  the  gage,  where  is  the 

molecule  mean  speed  at  the  temperature  of  the  gas  and  /YVj  is 
the  incremental  molecule  number  density.  Of  these 
succeed  in  emerging  from  the  entrance  tube,  where  W3  is  the 
probability  that  a molecule  striking  the  area  05  b will 
emerge  through  the  area  Gy),  without  re-entering  the  gage. 

Wj  depends  upon  the  shape  of  the  entrance  tube  only. 

At  equilibrium,  then. 


<2-3U) 

Since  this  adjustment  takes  place  under  free  molecule 
flow  conditions  it  is  assumed  to  be  independent  of  other  simi- 
lar adjustments  taking  place  simultaneously  as  a result  of 
small  changes  in  instrument  base  pressure,  for  example.  It  has 
also  been  assumed  that  no  thermal  gradients  exist  along  the 
entrance  tube. 


In  order  to  fix  ideas  of  the  orders  of  magnitude  involved 
in  the  detection  of  a diffusely  scattered  beam,  consider  the 
idealized  schematic  of  a scattering  experiment.  Figure  7. 

When  the  length  of  the  beam  path,  A , very  greatly  ex- 
ceeds cl,  , the  radius  of  the  source  orifice,  point  source 
geometry  affords  a good  approximation  in  the  calculation  of 
beam  intensity.  The  number  of  molecules  effusing  from  the 
source  orifice  under  conditions  of  free  molecule  flow  is 
/n^»  0*  , where  0r  is  the  area  of  the  source  orifice, 

C 1 is  the  mean  speed  of  the  source  molecules,  and  rt\,  is 
the  source  molecule  number  density.  Let  the  beam  cross  sectional 
area  be  chosen  as  Ciz  at  A and  the  defining  slit  size  adjusted 
to  permit  this.  The  number  of  molecules  crossing  0*  per  second 
becomes 


No  = mi  c,a»(u 

8 tta*- 


(2.312) 
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Let  these  molecules  be  intercepted  hy  the  scattering  surface. 

For  an  unanbiguous  determination  of  the  flux  of  molecules  through 
the  area  Q3  at  angle  0'  as  measured  from  the  normal  to  the 
scattering  surface,  let  Cl 3 be  an  orifice  in  the  mall  of  the 
detector  gage  ( W*  1)  . 

If  there  is  perfectly  diffuse  scattering  at  the  surface, 
the  number  of  molecules  crossing  O3  into  the  gage  is 
£j_  & t S & # At  equilibria  then 

nr>  3 q-x  = c,  fl.  (X-L  Cx-t  coi  Q>x  since  -Sp-a  YS 

A 16  Trv>'iV-  Ci  'T> 

and  this  ratio  is  of  order  1 for  moderate  gage  temperatures. 

Thus: 

cos  Q*  . . 

+/?  R*  12.313) 

where  rv^is  the  incremental  molecule  number  density  in  the  gage 
due  to  the  scattered  beam. 

The  desired  angular  resolution  of  the  experiment,  the  physi- 
cal size  of  various  components  of  the  beam  foiming  and  detect- 
ing system,  and  the  available  pumping  speeds  on  source  and  de- 
tector sides of  the  plane  of  the  defining  orifice  all  play  a part 
in  limiting  the  choice  of  the  parameters  on  the  right  hand  side 
of  equation  2.313.  For  the  present  experiment,  typical  values 

0.25  mm 
1.5  mm 
A * 25  cm 
* 3 cm 

11  -3 

^1=  3.6  x 10  molecules  cm  at  a maximum 
7 • 

then  ^3*  2.25  x 10  cos0  • The  pressure  increment  cor- 
responding to  this  density  increment  is 

Ap  — 0*0  3 xlO  /m/vn 


An  entrance  cone  of  taper  such  that  no  molecule  can  hit  the 
side  wall  before  entering  the  gage  volume  (Figure  7)  can  be 
constructed  to  have  a value  of  W3  of  the  order  0.2,  The 

scattered  beta  pressure  increment  may  thus  be  raised  to  the 

_q 

order  of  3.1  x 10  mm  Hg. 
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A pressure  increment  of  this  small  size  would  not  be 
difficult  to  detect  were  it  not  for  the  presence  of  the  back- 
ground gas.  This  residual  gas,  existing  in  the  detector  region 
at  a mean  pressure  of  10”^  mm  Hg,  suffers  rapid  random  fluctu- 
ations in  pressure  as  a result  of  uneven  diffusion  pump  operation* 
These  pressure  fluctuations  under  operating  conditions  are  of 
the  order  of  3 x 1C“®  mm  Hg.  This  gas  also  suffers  slower 
■drifts*  in  pressure  in  response  to  minor  changes  in  diffusion 
pump  efficiency  or  system  outgassing  rate.  These  changes  in 
turn  may  be  caused  by  changes  in  pump  heater  current,  ambient 
temperature,  cooling  water  temperature  or  cold-trap  temperature. 

A background  pressure  change  of  5 x 10  mm  Hg  may  occur  in  a 

-7  -7 

few  seconds,  of  10  mm  Hg  in  several  minutes  or  5 x 10  mm  Hg 

in  the  course  of  a few  hours. 


Thus  the  atmospheric  constituent  detection  system  must  in- 
clude a second  gage  known  as  the  "comparator".  This  gage  must 
measure  the  background  pressure,  responding  to  its  fluctuations 
with  the  same  sensitivity  and  with  the  same  time  of  response  as 
that  of  the  detector.  If  the  comparator  provided  a perfect 
replica  of  the  signal  from  the  detector  in  response  to  the  back- 
ground pressure  fluctuations,  a subtraction  of  the  total  signal 
from  each  gage  would  leave  a net  signal  representing  only  those 
pressure  increments  due  to  the  beam.  This  condition  of  perfect 
compensation  is  never  quite  realized  in  practice. 


An  Important  requirement  for  good  compensation  is  that  the 
time  response  of  comparator  and  detector  gages  be  equal.  The 
time  response  of  either  gage  to  a stepwise  increment  of  mole- 
cular flow  into  the  gage,  M molecules  per  unit  time,  is  given 


by 


rv\(t) 


= JlU (i  -e 

C5W3&V 


g»w3q,t 

4 v ) (2.314) 


where  rv\(t)  is  the  incremental  number  density  of  molecules  in 
the  gage  volume  at  time  t after  initiation  of  the  incremental 
flow,  is  the  mean  molecular  speed  of  the  molecules  within 
the  gage,  Q3  is  the  area  of  the  entrance  tube  opening,  and 
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I,  is  the  probability  that  a molecule  will  pass  on  out  of  the 
tube  without  returning  to  the  gage  volume.  V is  the  volume  of 
the  gage.  If  N is  due  to  a stepwise  increase  in  molecule  number 
density  nrt o within  the  vacuum  shell  of  the  detector  region, 

(t)  =.  <v\0  C~»  ( i - 6.  * v V 

£3  ' ' (2.315) 


Clearly,  if  the  response  times  of  detector  and  comparator 
are  to  be  made  equal,  the  factors  ^ W3  (X-$  must  be 
made  equal.  This  was  most  conveniently  achieved  by  providing 
the  comparator  gage  with  a telescoping  entrance  tube  which 
was  then  adjusted  by  trial  and  error  for  minimum  time  response 
difference. 

Although  the  fundamental  concepts  of  the  detector-com- 
parator form  of  measurement  have  been  widely  applied  in  other 
fields  (and  referred  to  in  other  terms)  they  are  believed 
to  have  been  first  applied  to  a molecular  beam  detection  system 
by  0.  Stem,  I.  Estermann  and  their  associates. 

2.32  The  Ionization  Gage  Detector 

As  has  been  remarked  earlier  in  this  report,  the  ioni- 
zation gage  was  chosen  for  development  as  a molecular  beam 
detector.  The  Plrani  gage  appeared  to  have  been  brought  close 
to  the  limit  of  its  potential  by  earlier  workers  (again  Ester- 
mann, Stern,  et.  al).  It  did  not  seem  that  their  work  could 
be  bettered  by  one  or  two  orders  of  magnitude. 

In  its  most  common  form  the  ionization  gage  employs  three 
major  elements,  contained  within  and  supported  by  a glass  envelope. 
The  heated  filament,  operating  temperature -limited,  supplies 
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electrons  for  the  ionisation  of  neutral  molecules  in  the  grid- 
anode  region*  The  grid  accelerates  the  electrons  and  ultimately 
collects  them  after  they  have  made  several  passes  from  grid- 
cathode  to  grid-anode  region  and  hack*  Such  ions  as  are  formed 
hy  electron  collision  in  the  grid-anode  region  are  collected 
at  the  anode*  Electrode  potentials  and  field  configuration 
are  such  as  to  permit  this  action*  The  magnitude  of  the  ion 
current  is  linearly  proportional  to  the  number  density  of 
molecules  in  the  gage  at  pressures  below  one  micron  Hg  and  to 
the  electron  emission  current.  It  is  also  a function  of  grid 
potential  and  gage  geometry,  and  is  inherently  very  insensitive 
to  changes  in  ambient  temperature. 

An  emitting  tungsten  filament  continuously  loses  mass  to 
its  surroundings*  Ultimately,  conducting  films  of  tungsten 
be  cane  deposited  on  the  gage  walls*  At  the  same  time  diffusion 
pump  oils  are  broken  down  to  elemental  constituents  or  to 
resinous  molecules  and  also  become  deposited  on  the  walls* 

A complex  film  is  built  up  having  conduction  properties  of 
varied  and  sometimes  rapidly  fluctuating  character*  Leakage 
paths  and  fluctuating  charge  distributions  have  a very  dis- 
turbing effect  on  the  stability  of  the  ion  current  when  that 
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current  must  be  known  to  a few  parts  in  10  , although  for  or- 
dinary pressure  monitoring  these  deposits  are  of  little  interest* 

In  designing  an  ion  gage  for  use  as  a molecular  beam  de- 
tector, particular  care  must  be  taken  to  insure  that  electrical 
leakage  paths  between  the  various  elements  or  from  them  to 
ground  have  virtually  no  opportunity  of  formation.  The  gage 
must  also  be  exceptionally  rugged  and  physically  compact 0 

By  February  of  1949,  the  major  design  features  of  a satis- 
factory ion  gage  for  molecular  beam  detection  were  well  estab- 
lished. A detailed  record  of  this  development  is  presented  in 
Ref*  26*  The  ion  gage  had  improved  itself  in  two  experimental 
programs  (Refs*  27  and  28)  and  it  had  bear  established  that 
signals  at  least  as  small  as  those  encountered  in  the  work  of 
earlier  experimenters  could  be  resolved  from  the  noise* 
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In  1952  a major  change  mas  made  in  the  molecular  beam 
detector  as  a result  of  the  publication  of  the  work  of  D.  Alpert 
and  others  at  Westinghouse  Research  Laboratory  (Ref.  29).  While 
the  major  advantage  of  the  Alpert  gage  (i.e.  the  measurement 
of  extremely  low  absolute  pressures)  was  thought  not  to  be  of 
interest  at  this  time  its  configuration  allowed  a cleanness 
and  ruggedness  of  structure  which  was  regarded  as  very  desirable* 

On  the  basis  of  the  Alpert  design  and  on  the  previously  achieved 
design  of  this  project,  the  ion  gage  detector  of  the  present 
.investigation  was  constructed  (see  Figure  8). 

The  sensitivity  and  other  performance  characteristics  of 
these  gages  were  discovered  to  be  very  similar  to  those  of  the 
Alpert  gage.  In  spite  of  the  attention  given  to  the  prevention 
of  film  formation  this  still  constituted  a serious  problem* 

One  month  was  found  to  be  about  the  maximum  running  time  be- 
fore conducting  films  made  the  gage  inoperable  for  molecular 
beam  use.  A thorough  cleaning  with  dilute  hydrofluoric  acid 
was  found  to  be  necessary  to  restore  a dirty  gage  to  useful- 
ness. 

2.33  ifag  mggjgB  JteoflateE 

For  the  proper  operation  of  the  ionization  gage  detector 
a stable  source  of  gird  potential  must  be  supplied.  Also  a 
source  of  filament  power  must  be  provided  which  is  so  regulated 
as  to  permit  an  unvarying  flow  of  emission  current  within  the 
gage.  It  was  a requirement  of  the  experiment  that  this  current 
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be  regulated  to  a few  parts  in  10  * Similarly  the  source  of 
grid  potential  was  required  to  possess  a very  high  order  of 
stability*  The  device  which  supplied  both  filament  and  grid  was 
termed  the  emission  regulator.  Two  such  devices  were  constructed, 
one  for  the  detector  gage  and  one  for  the  comparator.  (Beam 
detector  system  schematic.  Figure  9)* 

This  filament  emission  regulator  was  of  a type  in  which 
the  total  emission  current  develops  a signal  «cross  the 
control  resistor  R*  (Figure  10)*  The  potential  E^  is  compered 
with  a reference  potential  EQ.  The  difference  eQ  (E1-Eq  ■ e^) 
constitutes  an  error  signal  whish  is  then  minimized  by  the 
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regulator  action.  In  a total  feed-back  circuit  of  this  type 

A where  A is  the  loop  gain  of  the  circuit.  All  alnor 
Eq  1 + A 

fluctuations  in  regulator  performance,  supply  voltage,  filament 
emissivity, etc.,  can  be  considered  equivalent  to  changes  in  loop 
gain.  In  principle  the  ratio  jrj  can  be  made  to  remain  as  con- 
stant as  desired  by  providing  a sufficiently  high  loop  gain. 

In  particular  if  variations  in  effective  loop  gain  of  10$  A are 

Ei 

to  be  expected,  and  if  the  ratio  g*  is  to  remain  constant  to 
1 part  in  10^  then  A must  have  the  value  of  10,000.  Vari- 
ations in  the  reference  potential  Eq  and  variations  in  input 
stage  quiescent  point  are  excluded  from  this  argument,  of  course. 

It  is  a matter  of  experimental  record,  however,  that  the 
achievement  of  a loop  gain  of  this  order  is  not  feasible  in 
the  straight  forward  D.C.  amplifier  type  of  emission  regulator 
design.  If  the  voltage  output  of  the  regulator  is  sinusoidal 
the  resultant  emission  current  will  have  a ripple  component 
superimposed  on  the  D.C.  level.  This  ripple  component  is  the 
result  of  the  response  of  the  temperature  limited  emission 
current  to  the  periodic  fluctuation  of  the  impressed  power. 

A high  gain  D.C.  amplifier  normally  becomes  blocked  by  this 
signal  in  a final  stage.  Attempts  to  reduce  the  magnitude  of 
this  ripple  component  by  brute  force  filtering  are  met  with 
very  limited  success  because  of  the  attendant  phase  shift  around 
the  loop  and  the  resultant  low  frequency  hunting. 

A point  at  which  phase  shift  is  unavoidably  introduced  is 
at  the  filament  of  the  ionization  gage  itself  as  a result  of  the 
thermal  capacity  of  the  filament.  This  highly  non-linear  ele- 
ment makes  the  realization  of  extremely  precise  regulation 
more  difficult  than  in  the  case  of  a pure  resistive  load. 

The  design  of  the  ion  gage  detector  regulator  followed 
an  approach  suggested  by  Greenough,  Williams  and  Taylor  (Ref.  30) 
and  applied  by  them  to  a low  voltage  power  supply  for  electro- 
lysis. The  essential  feature  of  the  Greenough,  Williams  ap- 
proach was  the  use  of  a square  wave  as  the  regulator  output 
signal.  The  application  of  power  in  this  form  to  the  filament 
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result*  in  a great  reduction  in  the  eniesion  current  ripple 
oonponent  and  a consequent  increase  in  the  feasible  loop  gain. 

Schematic  diagrams  of  the  complete  regulator  and  power 
supply  circuits  are  presented  in  figures  11  and  12.  Referring 
to  the  schematic  of  Figure  11  the  emission  current  entering 
at  point  A develops  a signal  voltage  across  precision  resistors 
R.  This  voltage  is  compared  with  reference  battery  EQ  and  the 
resultant  difference  signal  is  applied  to  the  grid  of  the 
first  amplifier  Vg.  compensates  the  signal  for  small  changes 
in  contact  potential  originating  from  filament  emission  changes 
in  Vg.  V4  and  V,.  comprise  a multivibrator  pair  serving  as  the 
plate  load  of  the  control  tube  ly  A square  wave  of  acceptable 
quality  is  electron  coupled  out  through  a loop  gain  control 
dual  potentiometer  to  a high  grain  AC  amplifier  adapted  from 
the  Childs  circuit  (Ref.  31).  The  output  transformer  T^  was 
chosen  to  have  db  band  pass  of  20-20,000  cycles,  to  Insure  good 
transfer  of  the  square  wave  to  the  ion  gage  filament.  Provision 
was  made  for  balancing  the  plate  currents  in  the  push-pull 
power  output  stage.  It  was  found  necessary  to  include  the 
variable  RC  network  R^  R2  R^  C^  Cg  as  an  anti -hunt  control. 
Highest  performance  of  the  regulator  was  achieved  by  adjusting 
anti-hunt,  loop  gain,  and  various  balance  controls. 

The  power  supply  units  (Figure  12)  were  rebuilt  from  Signal 
Corps.  SCR  545A  mobile  Radar  "standard11  power  supplies. 

Using  a precision  voltage  divider,  biasing  batteries  and 
a 0-20  millivolt  strip  chart  recorder,  determinations  of  the 
stability  of  the  regulator  and  of  its  loop  gain  were  made. 

After  a short  initial  warm  up  period  the  "drift"  was  of  the  order 
of  one  part  in  104  per  minute.  This  was  a matter  of  no  concern 
since  after  warm  up  the  drifts  had  been  ascertained  to  be  less 
than  1%  per  day  except  under  conditions  of  excessive  ambient 
temperature  change.  (The  reference  potentials  were  supplied 
by  ordinary  dry  cells) • The  noise  was  of  the  order  of  2 parts 
in -10'’  of  the  total  signal  in  amplitude,  permitting  adequate 
gage  operation.  The  loop  gain  was  variously  checked  at  between 
10,000  and  20,000  under  good  operating  conditions. 
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While  long  tine  stability  of  the  highest  order  ns  not  a 
requirement  of  the  program  of  this  report,  it  could  have  been 
approached  more  closely  by  making  two  simple  changes  in  the 
equipment.  The  reference  potentical  could  have  been  supplied 
by  a mercury  cell  having  a low  temperature  coefficient  of  e.m.f., 
and  the  heaters  of  the  tubes  of  the  input  stage  could  have  been 
powered  from  the  regulated  supply. 

2.34  Iha 

At  some  point  in  the  detection  process,  in  order  that 
meaningful  data  can  be  obtained  from  the  molecular  beam  detector, 
the  signals  from  detector  and  comparator  must  be  compared  and  a 
difference  signal  noted.  The  apparatus  which  performed  this 
function  in  the  molecular  beam  detector  system  was  known  as  the 
difference  amplifier. 

In  order  to  understand  the  nature  of  the  problem  it  is 

useful  to  recall  the  orders  of  magnitude  of  the  current  involved. 

—8  “7 

The  background  pressure  current  was  of  the  order  of  10  to  10 

„9 

amperes.  Random  current  fluctuations  of  the  order  of  10  amperes 
in  amplitude  wars  superimposed  on  this  as  a result  of  itoe  background 
pressure  fluctuations.  The  beam  signal  current  was  of  the  order 
of  10-10  amperes.  It  was  thought  desirable  to  resolve  this  sig- 
nal to  at  least  one  part  in  10.  The  conditions  of  the  experi- 
ment required  therefore  that  currents  from  detector  and  com- 
parator be  amplified  and  subtracted  in  such  fashion  as  to  yield 
difference  signals  of  as  little  as  one  part  in  10^  of  the  total. 

Systems  for  performing  a similar  function  have  been  devised 
in  the  past  for  use  in  mass  spectrometer  Isotope  abundance  stu- 
dies (see  for  instance.  Ref.  32).  These  systems  normally  used 
one  oomplete  amplifier  for  each  signal,  the  subtraction  being 
performed  directly  proceeding  the  indicating  device.  It  was 
thought,  in  view  of  the  special  circumstance  of  having  nearly 
identical  input  signals,  that  the  two  amplifiers  could  be  joined 
in  a single  push-pull  circuit  minimising  the  opportunity  for 
unequal  phase  shift  through  two  amplifiers  and  for  other  inequali- 
ties of  response  possible  with  two  oomplete  amplifiers.  This 
and  more  specialised  design  considerations  led  to  the  difference 
amplifier  circuit  shown  in  Figure  13.  Power  supplies  for  this 
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unit  are  shown  in  Figure  14. 


If  circuit  sections  of  Figure  13  including  V^,  V2  and 
V3,  are  considered  alone,  it  can  be  seen  that  these  together 
form  a basic  feed  back  electrometer  circuit  for  which  eQU^.  ■ 
i^n  R^  ( 2+^  ) eQu^  is  the  signal  potential  to  the  meter- 

ing circuit,  1^  is  the  signal  current  in,  is  the  high  re- 
sistance range  determining  resistor  and  A is  the  effective  gain 
of  the  amplifier  section.  As  with  the  emission  regulator,  if 
gain  variations  of  10$  A are  to  be  expected,  A must  have  a 


value  of  the  order  of  10^  for  the  ratio 


out 

"in*! 


to  remain  con- 
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stant  to  one  part  in  10  * If  a shunting  resistance  to  ground  is 
present  (as  for  instance  in  Figure  14)  the  product  A R2  must 
be  very  much  greater  than  R^  in  order  for  this  shunt  to  have  a 
negligible  effect  on  the  output  signal.  The  gain  of  the  ampli- 
fier in  the  difference  amplifier  circuit  was  somewhat  greater  than 
104. 


Referring  again  to  the  schematic  of  Figure  13,  the  signals 

from  detector  and  comparator  are  conducted  to  V_  and  V_  housed 
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in  the  pick  up  head.  This  section  of  the  unit  wa3  adapted  with 
minor  modification  from  the  National  Research  Corporation  Alpha- 
tron  vacuum  gage.  Selected  954  type  vacuum  tubes  connected 
as  shown  have  been  found  to  show  grid  currents  of  less  than  10”*4 
amperes  (Ref.  33).  The  pairs  V2,  and  V^,  were  connected 
with  common  cathode  resistors  to  minimize  the  effects  of  variation 
in  plate  supply  voltage  and  to  afford  a means  of  balancing  the 
gain  of  each  stage.  The  balanced  cathode  follower  output  stage 
provided  a low  impedance  drive  for  the  strip  chart  recorder. 

8 

With  an  input  range  resistor  of  10  ohms  the  least  output 

difference  signal  distinguishable  on  the  strip  chart  recorder 

-12 

was  of  the  order  of  10  amperes.  Drift  and  noise  appeared  to 
be  absent  at  this  sensitivity.  Larger  values  of  the  input 
range  resistor  were  not  employed  because  of  the  magnitude  of 
the  voltage  which  would  have  been  developed  across  than  by  the 
large  background  pressure  current* 

All  heaters  were  supplied  from  the  regulated  supply  of 
Figure  14.  This  supply  was  also  modified  from  the  Signal  Corps. 
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SCR545A  "standard*  power  supply. 

3.0  THE  EXPERIMENTAL  PROGRAM 

3.1  Test  Surfaces  and  Geometries 

In  view  of  the  survey  nature  of  this  experimental  program,  it 
was  regarded  as  of  first  importance  to  choose  test  surface  and  beam 
geometries  in  such  a way  as  to  allow  an  overall  determination  of  the 
scattering  behavior  with  a minimum  of  ambiguity.  The  results  of  the 
first  design  choices  are  shown  in  Figure  15.  The  beam  was  axially 
symmetric  as  was  the  gage  entrance  cone.  The  surface  was  circular 
and  of  approximately  the  same  area  as  that  of  the  entrance  end  of 
the  gage  entrance  cone.  At  normal  incidence  the  beam  umbra  illumi- 
nated somewhat  less  than  the  total  surface  area.  However,  low  inten- 
sity elements  of  the  penumbra  overlapped  the  surface  slightly.  The 
gage  entrance  cone  was  designed  to  accept  approximately  1%  of  the  total 
hemisphere  solid  angle.  This  was  a "low  resolution*  geometry,  chosen, 
on  the  basis  of  beam  intensity  calculations  and  available  signal-to- 
noise  ratio,  so  as  to  provide  a reasonable  expectation  of  yielding 
detectable  signals  for  the  case  of  diffuse  scattering. 

Three  surfaces  were  tested  using  this  test  geometry.  These 
surfaces  and  the  test  surface  holder  are  illustrated  in  Figure  16.  A 
summary  of  this  portion  of  the  testing  program  follows. 

3.11  tolled,,  S.tg.gJ, 

This  specimen  was  of  common  cold  rolled  steel  of  low  car- 
bon content,  lathe  turned  from  a l/4"  diameter  rod.  The  speci- 
men was  mounted  in  lucite  and  polished  using  standard  metal- 
lurgical laboratory  techniques.  Examination  revealed  a uniform 
mat  of  fine  polish  marks  of  the  order  of  l/4  micron  in  depth. 

The  surface  appeared  to  the  naked  eye  to  have  a high  polish. 

After  insertion  into  the  specimen  holder  in  the  molecular  beam 
apparatus,  the  surface  was  thoroughly  washed  with  absolute  al- 
cohol and  dried  with  tank  nitrogen.  Polar  traverses  were  made 
both  in  the  incident  plane  and  out  of  the  incident  plane  with 
the  beam  incident  at  a variety  of  representative  incident  angles. 
Room  air  served  as  the  test  gas.  The  surface  temperature  was 
very  slightly  above  room  temperature.  About  250  data  points 
were  obtained  in  the  period  Feb.  12,  1953  to  March  2,  1953. 
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Following  the  completion  of  this  sequence  the  surface  w 
heated  to  250  degrees  centigrade  and  a series  of  incident  pi 
polar  traverses  were  made  at  this  temperature  in  order  to  di 
cover  any  possible  change  in  the  scattering  behavior  as  a re 
of  the  outgasslng  of  organic  films  or  otiier  adsorbed  layers. 

A final  polar  traverse  was  made  again  at  room  temperatu 

3.12  Gold  Rolled  Steel.  Etched 

A specimen  of  the  same  material  as  that  of  item  3.11  wa 
polished  in  the  same  fashion.  It  was  then  given  a light  etc 
in  a standard  10$  nitol  solution.  Microscopic  examination  s 
the  characteristic  grain  structure  of  a low  carbon  steel.  N 
polish  marks  were  visible.  Room  air  was  again  used  as  the  b 
gas.  This  specimen  was  tested  as  before,  first  at  room  temp 
erature,  then  at  elevated  temperature,  and  again  at  room  tem 
ture.  Incident  plane  polar  traverses  only  were  made.  It  wa 
thought  not  worthwhile  to  spend  the  time  on  the  exhaustive  t 
of  the  first  section  unless  significant  changes  in  the  scatt 
ing  behavior  developed.  In  all,  about  200  data  points  were 
tained  in  the  period  March  5,  1953  to  March  13,  1953. 

3.13  Almlm 

This  specimen  was  lathe  turned  from  S02,  a commercially 
grade  of  aluminum.  It  was  polished  in  the  same  manner  as  th 
steel  specimen.  This  specimen  was  tested  as  in  item  3.12. 
together,  some  120  data  points  were  obtained  in  the  period 
March  23,  1953  through  April  8,  1953. 

Examination  of  the  data  of  these  first  tests  revealed  t 
deficiences  of  the  experimental  geometry  which  were  corrects 
in  part  before  proceeding  with  the  glass  surfaces.  It  was  o 
served  that  at  angles  of  incidence  moderately  close  to  glanc 
the  edge  of  the  specimen  intercepted  a portion  of  the  beam, 
detector  positions  back  of  the  test  surface  normal,  the  edge 
scattered  molecules  contributed  in  part  to  the  detector  sign 
making  difficult  an  interpretation  of  the  scattering  data  in 
those  directions.  Secondly,  direct  beam  molecules  were  dete 
at  incident  and  exit  angles  both  close  to  grazing  as  a resul 
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of  finite  beta  breadth*  While  such  conditions  are  ultimately 
met  in  any  beam  geometry  it  mas  thought  that  in  this  geometry 
the  limitation  on  large  incident  angle  measurement  was  excessive. 

The  correction  of  both  of  these  imperfections  could  have 
been  accomplished  by  changing  to  a narrow  ribbon  geometry  but 
only  at  the  sacrifice  of  out-of-the-incident-plane  information. 

It  was  thought  better  to  retain  the  overall  survey  approach  at 
this  stage,  but  to  narrow  the  beam  by  replacing  the  circular 
secondary  defining  orifice  by  a 1 mm  x 5 mm  slit.  It  had  been 
ascertained  that  the  signal  to  noise  ratio  was  sufficiently 
high  to  permit  this  attenuation  of  the  beam.  At  the  same  time, 
edge  scattering  was  minimized  by  increasing  the  test  surface 
length.  This  change  also  made  the  preparation  and  mounting  of 
the  glass  surfaces  less  difficult.  This  beam  geometry  is  illus- 
trated in  Figure  15.  The  glass  test  surfaces  are  included  in 
Figure  16. 

3.14  Unpolished  Glass  Surface 

This  specimen  was  obtained  from  a pane  of  single  strength 
window  glass.  It  was  thoroughly  cleaned,  given  a final  rinse  in 
absolute  alcohol,  and  dried  with  tank  nitrogen.  The  surface  as 
viewed  by  the  naked  eye  showed  no  defects  and  appeared  free  of 
film.  The  structure  of  the  surface  was  presumably  "as  cooled". 

It  was  this  specimen  which  was  chosen  for  the  second  ex- 
haustive scattering  examination,  the  results  of  which  were  later 
to  be  used  in  the  determination  of  f,  the  coefficient  of  specular 
reflection.  This  choice  was  in  part  determined  by  the  reported 
value  f * 0.89  for  air  on  glass  as  remarked  in  Section  1 of 
this  report.  In  all,  during  the  period  April  17,  1953  to  June 
21,  1953,  more  than  400  data  points  were  obtained.  Tank  nitro- 
gen (99.5%  Ng)  mas  used  as  the  test  gas.  The  tests  were  run 
with  the  surface  essentially  at  ambient  temperature.  Since  it 
was  mounted  on  a rod  which  was  warned  by  conduction  through 
parts  connected  to  the  detector  ion  gage,  the  test  surface  was 
maintained  at  about  10#C  above  the  temperature  of  the  surrounding 
walls. 
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3.15  M&t&qQlflgg  sarfaffift 

This  specimen  was  cut  from  the  same  stock  as  the  unpolished 
glass  surface.  It  was  then  sent  to  a manufacturing  optical  lab- 
oratory for  polishing.  The  polished  specimen  was  cleaned  and 
rinsed  In  the  same  fashion  as  before.  For  purposes  of  comparison 
with  the  unpolished  glass  plate  it  was  necessary  to  observe  scat- 
tering in  the  incident  plane  only.  Approximately  100  data 
points  were  obtained.  The  experimental  program  of  this  report 
was  completed  with  the  conclusion  of  these  tests  July  23,  1953. 

3.2  The  Experimental  Data 

Each  polar  traverse  made  at  fixed  incident  and  azimuth  angles 
yielded  several  data  points  representing  the  flux  of  scattered  mole- 
cules into  the  detector  gage.  Each  datum  point  resulted  from  the 
averaging  of  a number  of  beam  on-off  traces  as  recorded  by  the  strip 
chart  recorder.  A reproduction  of  such  a record  is  presented  in 
Figure  17.  For  each  day's  run  the  beam  source  pressure  was  held  con- 
stant as  was  the  gage  emission  current,  and  hence  the  overall  system 
sensitivity.  When  possible,  in  each  run  a record  of  the  normal  in- 
cidence scattering  was  taken  and  also  a record  of  the  direct  beam 
strength.  These  last  records  afforded  a means  of  comparison  of  one 
run  with  another. 

It  was  found  most  convenient  to  represent  the  data  in  polar  plot 
form,  scaling  the  data  of  each  run  up  or  down  to  bring  data  of  all 
runs  of  a given  sequence  into  register. 

In  the  case  of  the  steel  and  aluminum  surfaces  two  additional 
problems  in  data  interpretation  must  be  mentioned.  As  a result  of 
the  width  of  the  shutter  mechanism  and  of  the  detector  ion  gage  the 
polar  traverses  could  be  continued  to  within  about  40*  of  the  reverse 
beam  axis  only.  This  conditions  was  improved  only  slightly  in  the 
glass  plate  tests.  In  these  it  was  possible  for  the  detector  axis 
to  came  within  30  degrees  of  the  reverse  beam  axis.  The  second  pro- 
blem resulted  as  a consequence  of  the  small  test  surface  size  in  re- 
lation to  the  beam  diameter.  As  the  angle  of  incidence  was  increased, 
an  increasing  percentage  of  beam  molecules  failed  to  strike  the  sur- 
face. The  intensity  of  illumination  per  unit  of  surface  area  thus 
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fell  off  and  in  consequence  the  signal  amplitude* 

Zn  view  of  the  bulk  of  the  experimental  data  and  its  general 
similarity  in  the  case  of  the  steel  and  aluminum  surfaces  no  attempt 
will  be  made  to  present  it  in  its  entirety.  However,  representative 
polar  plots  for  the  steel  plate  are  given  in  Figures  18  through  21 
and  for  the  aluminum  plate,  Figures  22  through  25. 

Referring  to  the  figures  for  the  steel  surface  one  can  see  un- 
fortunate gaps  appearing  in  the  data  at  each  angle  of  incidence. 
However,  these  gaps  are  somewhat  filled  by  data  obtained  at  other 
angles  of  incidence.  The  most  striking  aspect  of  these  plots  is  the 
remarkable  adherence  of  the  data  to  a cosine  distribution.  The  plots 
of  data  taken  out  of  the  incident  plane  still  exhibit  the  same  degree 
of  circularity.  No  discernible  change  in  the  character  of  the  scat- 
tering was  made  by  outgassing  the  surface  or  running  traverses  at 
elevated  temperature.  For  these  surfaces  undistorted  representations 
of  forward  scattering  were  obtained  to  incident  angles  of  about  70 
degrees.  The  back  scattering  data  at  this  incident  angle  was,  how- 
ever, seriously  distorted  by  edge  scattering.  At  larger  angles  of 
incidence  (to  82  degrees)  no  evidence  appeared  of  marked  specular 
scattering,  although,  has  such  occurred  close  to  90  degrees,  detection 
would  have  been  masked  by  the  signal  of  the  direct  beam. 

A simple  comparison  of  these  data  with  the  reference  circles 
(introduced  to  aid  in  bringing  the  plots  into  register)  suggests  that 
the  amplitude  resolution  of  the  experiment  under  these  conditions  was 
of  the  order  of  two  to  three  percent  of  the  maximum  diffusely  scattered 
signal. 

Considering  now  Figures  22  through  25  for  the  aluminum  surface 
one  can  observe  that  most  of  the  remarks  applying  to  the  discussion  of 
steel  surfaces  might  also  be  applied  to  the  aluminum.  This  work  was 
not  as  complete  as  that  on  the  steel  surfaces,  being  done  in  the 
incident  plane  only  and  at  fewer  angles  of  incidence.  While  there 
seemed  to  be  some  deterioration  in  the  signal  amplitude  resolution 
over  this  period,  one  must  infer  that  the  data  indicates  only  slight 
deviations  if  any  from  the  cosine  distributions.  There  does,  however, 
seem  to  be  an  indication  of  small  bulges  in  the  plots  of  the  data  of 
March  23,  1953,  incident  angles  53.65*  and  81.83*. 
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The  plots  of  the  distribution  of  flux  from  the  glass  surfaces 
show  a somewhat  different  character#  Representative  data  are  illus- 
trated in  Figures  26  through  32  for  the  unpolished  glass  surface  and 
Figures  33  through  40  for  the  polished  glass  surface.  On  these  plots 
one  can  readily  observe  bulges  in  the  forward  scattering  direction, 
increasing  in  importance  as  the  angle  of  beam  Incidence  increases. 

These  bulges  are  of  the  same  general  magnitude  and  shape  for  the  two 
surfaces.  Thus  there  would  appear  to  be  some  interesting  but  un- 
drama  tic  deviations  from  the  pure  cosine  distribution.  While  scat- 
tering data  was  taken  at  an  angle  of  beam  incidence  as  great  as 
84.82*  no  large  specular  peaks  were  observed.  The  plots  at  this  angle 
of  Incidence  show  the  influence  of  edge  scattering  and  of  direct  beam 
detection  at  the  extremes  of  the  polar  traverse. 

One  further  obstacle  to  the  quantitative  interpretation  of  the 
data  appeared  in  the  observed  flattening  of  the  scattering  patterns 
at  the  poles.  This  flattening  arises  as  a consequence  of  the  elon- 
gation of  the  beam  spot  on  the  extended  test  surface  at  large  angles 
of  incidence.  Where  the  gage  entrance  tube  also  makes  a large  angle 
with  the  test  surface  normal  the  illuminated  spot  is  entirely  ■seen0 
by  the  gage.  As  the  tube  approaches  normal  incidence,  however,  the 
gage  can  only  "see*  a portion  of  the  illuminated  area.  If  completely 
diffuse  scattering  were  present  the  entire  pattern  would  lie  within 
the  scattering  circle  established  at  normal  incidence  and  would  be 
flattened  at  the  pole  but  symmetrically  disposed  about  the  test  sur- 
face nonnal. 

In  order  to  compare  the  phenomenological  results  of  the  scatter- 
ing tests  with  the  reported  value  of  f for  air  on  a glass  surface, 
an  attempt  was  made  to  organize  the  data  for  this  surface  into  com- 
plete "shells",  each  shell  representing  the  distribution  of  scattered 
flux  for  all  values  of  the  exit  coordinates  at  one  angle  of  incidence. 
This  work  resulted  in  the  construction  of  the  smoothed  figures  of 
Figures  29  through  32.  The  outermost  line  of  each  figure  represents 
the  flux  as  measured  in  the  incident  plane.  The  inner  lines  re- 
present the  flux  in  the  other  meridian  planes,  taken  at  10  degree 
intervals. 

Since  the  data  was  quite  incomplete  in  many  areas  it  was  necessary 
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to  fill  in  the  shells  on  the  basis  of  certain  reasonable  assumptions. 

It  was  assumed  that  the  shell  area  directly  intersected  by  the  incoming 
beam  had  no  sharp  cavity  or  protuberance.  It  was  further  assumed  that 
the  "back  scattering"  section  of  the  shell  followed  the  cosine  distribu- 
tion. This  assumption  was  supported  by  observed  data  as  far  as  it  went. 
This  assumption  is  further  supported  by  Zahl  in  his  observations  of 
the  scattering  of  Hg  from  NaCl  cleavage  planes.  (Ref.  35).  This 
means  that  the  surplus  of  molecules  forming  the  bulge  on  the  "for- 
ward scattering"  side  were  assumed  to  be  withdrawn  from  the  total 
number  of  molecules  striking  the  surface,  leaving  the  remainder  to 
be  scattered  in  totally  diffuse  fashion.  Finally,  it  was  assumed  that 
the  fall-off  of  the  non-diffuse  component  away  from  the  incident  plane 
proceeded  in  the  same  gradual  fashion  as  did  the  observed  fall-off  of 
this  component  about  the  specular  ray  axis.  This  assumption  was  again 
supported  by  such  data  as  was  obtained  out  of  the  incident  plane. 

Four  such  shells  were  completed  in  detail  for  the  incident  angles 
53.65  degrees,  66.63  degrees,  79.62  degrees,  and  84.82  degrees.  Data 
for  incident  angles  at  0.00  degrees,  29.10  degrees  and  40.67  degrees 
yielded  completely  spherical  shells  or,  in  the  40.67  degree  case,  a 
shell  so  nearly  spherical  as  to  make  uneconomical  the  rather  laborious 
numerical  computations  to  be  described  in  section  4.0. 

3.3  A Discussion  of  Errors 

The  usability  of  the  data  of  this  experimental  program  was  not 
dependent  upon  a knowledge  of  the  absolute  value  of  a pressure  or  a 
current,  but  rather  upon  the  stability  of  the  experimental  condition 
over  the  period  of  one  data  sequence.  A data  sequence  would  be  any 
group  of  related  polar  traverses  including  some  trace  suitable  for  re- 
lating this  sequence  to  another.  It  was  also  dependent  upon  the  short 
time  stability  of  the  signal  over  the  period  of  a beam  on-off  cycle. 

The  factors  influencing  the  long  time  stability  of  the  experi- 
mental conditions  were  the  source  pressure,  the  ion  gage  sensitivity, 
and  the  difference  amplifier  sensitivity.  Of  these,  only  the  source  pres- 
sure variation  produced  a measurable  change  in  the  ultimate  signal  ampli- 
tude over  the  period  of  a run.  The  source  pressure  responded  somewhat 
to  ambient  temperature  changes  as  a result  of  differential  expansion 
in  the  source  needle  valve.  In  the  operating  range  (400  to  800  microns 
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Hg)  the  slope  of  the  source  pressure -beam  flux  curve  varied  from 
0,3  to  0.1  chart  units  per  micron.  However,  only  on  rare 
occasions  were  the  fluctuations  within  the  period  of  a data  sequence 
important  enough  to  warrant  a correction.  A one  percent  variation  in 
ultimate  signal  amplitude  from  one  end  of  a run  to  the  other,  was 
•typical  and  was  regarded  as  of  no  significane  for  the  quantitative 
or  interpretive  purposes  of  this  program. 

Less  can  be  said  about  changes  in  Inherent  gage  sensitivity 
except  that  these  changes  (apparently  due  to  film  formation  on  the 
gage  walls)  occurred  slowly  and  did  not  significantly  influence  the 
self  consistency  of  the  data. 

The  ultimate  assurance  of  satisfactory  long  time  system  stability 
arose  from  the  comparison  of  direct  beam  signal  strength  at  the  begin- 
ning and  end  of  a data  sequence.  Such  checks  showed  the  direct  beam 
signal  strength  variation  to  be  usually  smaller  than  one  percent  of 
the  total  beam.  Occasionally,  the  variation  was  no  greater  than  this 
over  a period  of  several  days. 

A much  more  serious  problem  was  that  of  achieving  a satisfactory 
short  time  fluctuation  level.  This  had  been  discussed  in  section 
2.3  on  instrumentation.  However,  it  is  worth  considering  the  experi- 
mentally observed  signal-to-noise  ratio  in  relation  to  the  reliabiliiy 
of  the  data.  In  the  normal  observing  range  the  difference  amplifier 
sensitivity  was  lO"^  amperes  per  chart  unit.  One  chart  unit  for 
this  recorder  was  equivalent  to  0.1  millivolt  and  occupied  approxi- 
mately 0.05  inches  of  the  recorder  scale.  At  normal  gage  sensitivity 
a deflection  of  one  chart  unit  would  be  produced  by  a pressure  change 
in  one  ion  gage  of  the  order  10-10  mm  Hg.  On  this  scale  the  typical 
direct  beam  signal  was  of  the  order  of  900  chart  units  while  the 
diffusely  scattered  signal  produced  a plotted  amplitude  circle  of 
the  order  of  18  chart  units  in  diameter.  This  is  consistent  with 
the  acceptance  solid  angle  of  the  detector  entrance  tube. 

At  times  when  the  observing  conditions  were  exceptionally  good 
( for  instance  as  represented  in  Fig.  17)  repeated  points  frequently 
reproduced  themselves  to  0.1  chart  units  and  polar  traverses  showing 
cosine  scattering  behavior  appeared  to  have  an  internal  consistency 
allowing  a scatter  of  0.1  to  0.3  chart  units  from  the  reference  circle. 
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Root  mean  square  deviations  of  the  signal  trace  amplitude  under  these 
conditions  irere  of  the  order  of  0.2  chart  units.  Observing  conditions 
were  regarded  as  too  poor  for  the  taking  of  in terpre table  data  when  the 
r m s deviations  exceeded  2.0  chart  units.  The  bulk  of  the  data  was 
obtained  under  conditions  lying  in  between  these  extremes.  Prom  ob- 
servation of  the  polar  traverses  it  was  judged  that  over  90%  of  the 
experimental  points  lay  within  one  half  chart  unit  of  a faired  curve. 
Least  square  curve  fittings  were  not  attempted,  since  it  was  clear 
that  the  shapes  of  the  scattering  plots  were  well  defined.  Unfortu- 
nately, there  remained  those  whose  scattered  experimental  points  led 
to  curve  fitting  judgements  which  may  have  been  in  error  by  as  much 
as  10£. 

3.4  Physical  Interpretation  of  the  Data 

The  salient  geometrical  features  of  the  data  may  be  reviewed  as 
follows : 

(1)  The  most  characteristic  shape  of  all  the  scattering  flux 
shells  was  that  of  a sphere.  In  the  case  of  the  glass 
surfaces  the  shells  showed  a bulging  distortion  of  the 
forward  hemisphere.  The  spherical  shape  (cosine  distribu- 
tion) was  nearly  perfect  in  the  case  of  the  metal  sur- 
faces tested. 

(2)  No  shell  possessed  a "lobe"  or  dramatic  protuberance. 

There  were  no  cusps,  cavities  or  reversals  of  curvature, 
i.e.,  the  local  curvature  was  always  inward. 

(3)  The  back  scattered  hemisphere  appeared  to  be  spherical 
surfaces  without  distortions  in  so  far  as  could  be  de- 
termined. 

(4)  The  distortion  of  the  shells  in  the  case  of  the  glass 
surfaces  increased  with  increasing  angle  of  incidence. 

The  point  of  maximum  distortion  occurred  close  to  the 
intersection  point  of  the  specular  ray. 

(5)  The  scattering  behavior  for  the  unpolished  and  polished 
glass  surfaces  was  very  similar. 

(6)  There  was  no  discerns ble  effect  on  the  scattering  distri- 
butions for  the  metal  surfaces  as  a result  of  outgassing 
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the  surfaces  either  during  or  following  the  period  of 
elevated  surface  temperature. 

These  data  provide  a confirmation  in  some  degree  of  what  is 
alreacty  known  or  what  might  be  guessed  regarding  the  interaction  of 
air  or  nitrogen  with  the  surfaces  tested.  That  is,  it  is  shown  that 
by  far  the  greatest  fraction  of  air  molecules  incident  to  these  sur- 
faces scatter  from  them  according  to  a cosine  distribution.  Values 
of  f for  air  on  steel  and  aluminum  have  not  yet  appeared  in  the 
literature.  On  the  basis  of  these  measurements  f = 1.0.  If  all 
scattering  shells  in  a sequence  were  figures  of  revolution  about  the 
polar  axis  but  not  necessarily  spheres,  f would  also  turn  out  to  be 
1.  It  is  intuitively  satisfying  to  conclude  that  probably  most  other 
surfaces  for  which  f equals  1 would  yield  cosine  distribution  scatter- 
ing shells  in  test. 

Still,  little  understanding  is  gained  of  the  basic  mechanism  of 
the  scattering.  It  was  not  discovered,  for  example,  whether  the 
interactions  of  the  single  molecules  at  the  surface  were  more  elastic 
or  more  inelastic  in  nature.  A perfectly  rough  surface  might  well 
give  rise  to  the  cosine  scattering  distribution  even  if  every  indi- 
vidual molecule-surface  interaction  were  perfectly  elastic.  One 
would  suspect  from  published  values  of  the  energy  accomodation 
coefficient  that  this  is  not  the  case.  On  the  other  hand,  one  would 
also  suspect  from  these  published  values  of  &<.  that  the  interactions  are 
not  entirely  inelastic.  For  example,  for  air  on  polished  cast 
iron  has  values  of  0.87  to  0.88,  although  f for  polished  cold  rolled 
steel  has  the  value  1.0.  However,  if  multiple  impacts  with  the  sur- 
face irregularities  are  not  excluded,  the  individual  interactions 
might  be  quite  elastic  and  still  result  in  the  above  values  of  f 
and  . (Ref.  34).  Thus,  while  the  results  of  these  experiments 
give  issuing  flux  distributions,  it  must  again  be  emphasized  that 
they  tell  us  little  about  the  energy  exchange  at  the  surface  and 
hence  little  about  the  individual  mole cule-surf ace  interaction  mech- 
anism. 

One  important  uncertainty  in  the  physical  interpretation  of 
these  results  is  that  relating  to  the  adsorbed  films  which  might 
lie  on  the  test  surface.  In  the  case  of  the  metal  surfaces  the  re- 
moval of  superficial  surface  films  (if  present)  by  heating  had  no 
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effect  on  the  test  results  as  noted  In  Item  (6)  above.  Since  the 
metal  surface  scattering  distributions  differed  from  those  of  the 
glass  surface  one  might  imagine  that  some  essential  properties  of  the 
underlying  surface  were  reproduced  in  the  adsorbed  films.  However, 
the  extent  to  which  the  basic  properties  of  the  surface  were  obscured 
ty  the  surface  films  remained  largely  undetermined  as  did  the  com- 
position or  extent  of  the  surface  films  themselves.  Since  oil 
diffusion  pumps  were  employed,  it  is  quite  possible  that  oil  mono- 
layers  formed  on  all  test  surfaces  in  spite  of  efficient  cold  trap- 
ping. It  is  also  quite  possible  that  nitrogen  and  other  residual  gas 
molecules  were  adsorbed  on  the  surfaces.  The  point  of  view  which  was 
therefore  adopted  as  a result  of  the  proceeding  consideration  was  that 
the  effects  of  surface  contamination  on  the  angular  flux  distributions 
were  small  (as  suggested  ty  the  results  of  item  (6)  above)  but  that 
some  reservation  must  be  retained  in  predicting  behavior  for  surfaces 
of  absolute  cleanness  on  the  basis  of  the  present  tests.  A modifi- 
cation of  the  apparatus  which  would  make  future  experiment  less  am- 
biguous in  this  respect  would  be  the  substitution  of  mercury  vapor 
pumps  for  the  oil  diffusion  pumps.  Mercury  vapor  does  not  possess, 
as  strongly  as  oil  vapor,  the  tendency  to  stick  to  metal  or  glass 
surfaces. 

The  scattering  behavior  evidenced  by  the  glass  surfaces  suggests 
strongly  that  a not-negligible  minority  of  the  incident  molecules 
suffered  elastic  impacts  with  tilted  surface  elements,  thus  smearing 
out  the  “specular”  component.  Since  the  majority  of  molecules  might 
still  be  considered  to  have  undergone  a cosine  scattering  one  might 
imagine  that  here  and  there  the  local  field  contours  have  joined 
to  form  elastic  scattering  regions,  standing  like  islands  in  a 
jumbled  sea  of  unjoined  fields.  This  model  of  the  inelastic  inter- 
action for  the  majority  of  incident  molecules  with  a smeared  out 
specular  reflection  for  the  remainder  is  interestingly  close  to  the 
original  Maxwell  model.  This  particular  picture  will  be  referred  to 
as  the  "modified  Maxwell  hypothesis". 

This  picture  is  by  no  means  the  only  one  permissible  as  an  in- 
terpretation of  the  observed  scattering  distributions.  As  was  re- 
marked earlier  in  this  section,  many  hypotheses  concerning  the 
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energy  exchange  between  molecule  and  surface  could  be  made 
which  would  be  consistent  with  the  observed  scattered  flux  distribu- 
tions * One  such  is  that  all  molecules  are  completely  accommodated  in 
energy  to  that  of  the  surface  and  are  re-emitted  from  the  surface 
with  a speed  distribution  in  a given  direction  as  though  they  had 
originated  from  a gas  in  equilibrium  at  the  temperature  of  the  sur- 
face. 

The  tests  on  the  polished  glass  surface  suggest  that  the  ele- 
mentary scattering  elements  retain  their  basic  character  in  spite  of 
the  polishing  operation.  Some  substantiation  of  this  is  found  in  the 
observations  of  lens  makers  that  the  polishing  operation  results  in 
some  local  melting  and  flow  of  the  glass  surface. 

Certainly  there  is  a vast  amount  of  work  to  be  done  before  much 
more  tinderstanding  of  the  gas  surface  interaction  with  such  surfaces 
as  these  could  be  claimed.  However,  one  of  the  most  promising  avenues 
of  future  discovery  would  appear  to  be  in  the  use  of  a velocity  se- 
lector. Such  a device  could  certainly  tell  whether  or  not  the 
"specular"  component  was  accommodated  to  the  surface  temperature,  and 
at  the  same  time  of  course,  it  could  yield  data  leading  to  an  evalu- 
ation of  the  accommodation  coefficient. 

There  are  other  experiments  also  likely  to  yield  interesting 
results.  Hhat  are  the  effects  of  raising  or  lowering  the  mean  in- 
cidlent  momentum?  At  what  angles  of  incidence  on  complex  polished 
surfaces  do  sharp  reflection  (if  any)  occur?  What  can  be  deduced 
about  complex  surfaces  and  their  associated  gas  surface  interactions 
from  careful  experiments  on  pure  single  ciystals?  These  questions 
suggest  only  a few  of  the  problems.  It  is  hoped  that  time  and 
further  experiment  will  provide  an  understanding  that  is  still  lack- 
ing with  the  completion  of  this  report. 

4.0  A QUANTITATIVE  EXAMINATION  OF  THE  DATA 
4.1  Introduction 

It  is  of  interest  to  consider  the  application  of  the  experimental 
scattering  distribution  to  the  calculation  of  f , the  coefficient  of 
specular  reflection.  In  performing  this  analysis  the  surface  under 
consideration  will  be  regarded  as  bounding  a hypothetical  free  mole- 
cule flow  having  mass  velocity  directed  tangentially  to  the  surface. 
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Incident  and  Issuing  momentum  calculations  mill  be  made,  subject  to 
various  assumptions  enumerated  below.  The  issuing  momentum  calcu- 
lations will,  in  part,  be  dependent  upon  the  form  of  the  experimental 
scattering  distributions.  From  the  incident  and  Issuing  momentum 
calculations,  values  of  f will  be  obtained  in  accordance  with  the 
definition  of  that  coefficient. 

The  following  assumptions  will  be  applied  to  these  calculations) 

A.  For  the  molecule-surface  interaction  the  element  of  surface  is 
assumed  to-be  isotropic,  l.e.,  molecules  entering  at  @ and  (J>  exhibit 
a characteristic  scattering  behavior  independent  of  (J>  (see  Fig.  41)# 
This  assumption  does  not  apply  to  cleaved  surfaces  of  regular  single 
crystals.  It  would  seem  reasonable,  however,  to  apply  it  to  the 
present  case  of  a cooled  liquid  surface  or  to  the  case  of  "man-finished” 
surfaces  of  micro  crystalline  or  amorphous  materials. 

B»  It  is  assumed  in  the  absence  of  experimental  evidence  to  the  con- 
trary that  the  form  of  the  scattering  distribution  is  dependent  on 
gas,  surface,  and  incident  angle  only  and  is  independent  of  the  speed 
of  the  Incident  molecules  within  the  range  of  speeds  possessed  by 
molecules  of  gases  at  ordinaiy  temperatures.  This  implies  that  the 
form  of  the  distribution  is  independent  of  the  distribution  of  speeds 
in  the  incident  ray.  It  may  very  well  be  that  this  assumption  will 
be  proved  inaccurate  by  later  experiment.  It  is  Introduced  here 
simply  to  help  provide  a model  for  the  calculations. 

C.  It  is  assumed  for  this  calculation  that  the  incident  gas  is  so 
rarefied  that  the  motion  of  the  impinging  molecules  is  essentially 
unaltered  by  those  scattered  from  the  surface  in  question.  The  velo- 
city  distribution  of  the  impinging  molecules  will  be  taken  to  be  the 
llaxwell  distribution  with  superimposed  mass  velocity* 

D.  The  definition  of  f will  again  be  stated,  i.e.,  that 


where  (Vt  is  the  magnitude  of  the  net  tangential  momentum  brought  up 
to  the  surface  per  unit  time  by  all  incident  molecules  and  <5*-  is  the 
magnitude  of  the  tangential  momentum  imparted  by  the  surface  to  the 
issuing  molecules  per  unit  time. 

At  this  point  two  alternate  assumptions  relating  to  the  mode 
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of  energy  accommodation  at  the  surface  nrast  be  introduced.  Calcu- 
lations made  on  the  basis  of  each  of  these  assumptions  will  be  pre- 
sented. 

£-1  It  is  assumed  that  molecules  issuing  in  any  particular  direction 
from  the  scattering  surface  possess  the  distribution  in  speed  which 
they  would  have  had  had  they  issued  from  a body  of  gas  in  thermal 
equilibirum  at  the  surface.  It  follows  that  the  thermal  accommodation 
coefficient  has  the  value  1. 

E-2  Consider  again  the  form  of  the  scattering  distribution  of  sec- 
tion 3.0.  As  was  more  fully  discussed  in  that  section,  the  scatter- 
ing shells  (represented  in  Figures  29  through  32)  show  a hemispheric 
shape  as  the  "back  scattering"  section.  The  distortions  from  the  cosine 
scattering  sphere  take  the  form  of  bulges  protruding  from  the  "for- 
ward scattering"  hemisphere.  It  is  assumed  that  those  molecules  con- 
tributing to  the  bulging  deviation  from  the  cosine  scattering  sphere 
have  been  perfectly  elastically  scattered  and  so  possess  the  same 
distribution  in  speed  as  before  impact.  The  remainder,  i.e.,  all 
those  contributing  to  the  cosine  scattering  sphere  are  assumed  to 
possess  the  distribution  in  speed  given  by  assumption  (E-l).  Thus, 
for  the  purpose  of  this  calculation,  the  great  majority  of  molecules 
are  assumed  to  have  been  diffusely  scattered  with  accommodation  coef- 
ficient 1 while  the  remainder  are  assumed  to  have  been  specularly 
reflected  from  tilted  surface  elements.  This  is  a special  case  of 
the  modified  Maxwell  hypothesis  discussed  earlier  in  this  report. 

4.2  Calculation  of  f Under  Terms  of  Assumption  E-l 

First  consider  the  quantity  of  tangential  momentum  incident  on 
the  surface  element.  For  a representation  of  surface  and  coordinate 
system,  see  Figure  41.  Let  the  element  of  Surface  be  placed  at  the 
origin  of  coordinates.  Let  the  magnitude  of  the  mass  velocity  be 
U • Let  this  mass  velocity  be  directed  along  the  positive  X axis. 
Consider  those  molecules  so  directed  as  to  strike  the  surface  from 
below.  (This  choice  of  coordinates  is  dictated  by  convenience  since 
Cz  will  be  positive  for  those  molecules  striking  the  plate  and  C^ 
and  0^  will  be  positive  for  molecules  having  a line  of  flight  di- 
rected into  the  first  octant).  C , C , and  C are  the  x,  y,  and 

dc  y z 

a components  of  the  total  molecular  velocity  relative  to  the  fixed 
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coordinate  system. 

The  number  of  molecules  of  all  velocities  striking  unit  area  per 
unit  time  can  be  expressed  as 


Fy  4CxiCY  AG 


(4.2) 


where  is  the  Maxwell  distribution  with  superimposed  mass  velocity. 


Fys  e 


_ /3*  4 ■+  ^z  ^ 


(4.3) 


where  n is  the  molecule  number  density  of  the  gas. 


and  V 


m 


A = j&*.  , 

lr% 

is  the  most  probable  speed.  Each 


molecule  brings  up  to  the  surface  a component  of  the  tangential 
momentum  in  the  direction  of  the  x axis  having  the  value  nO^.  Thus 
the  total  x component  of  tangential  momentun  brought  up  to  a unit 
area  of  surface  per  unit  time  by  molecules  of  all  velocities  is  given 
by 


G»(  s.  j / C*  CxFy  ACx  JCydCz. 

J&  / oq/oo  (4.4) 


Integration  yields 

Gl  — AfVx  'VtC  U - fWN  <V\  Vvw  U 

4 ~ amr 


(4.5) 


The  calculation  of  Gy,  the  magnitude  of  the  net  tangential 
momentum  transferred  from  the  surface  to  the  issuing  molecules,  re- 
quires a little  more  work  since  in  some  fashion  the  direct  experimental 
data  must  be  applied.  Consider  the  expression 

Or  = /^  Nt(ct,ei  ,<(>;)?, 

SCifi  $i£'rp*$r  (4.6) 

X CrCOl4>|,S\NOr<)CCACt<iU»l<ioOr 
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Here  P(  (C i ,©  (.  ,4>i  ,C  y,0r  ,Cj>r  ")  is  a distribution  function  de- 
noting the  probability  per  unit  solid  angle  and  speed  range  that 
a molecule  with  incident  speed  and  direction  ©i.  , (|>c  will  be 
scattered  with  speed  in  direction  ©y»  , (J>r  . Here 

also  N [ lCc  ) ©t , 4>i ) represents  the  number  of  incident, 

molecules  per  unit  speed  range  and  solid  angle  having  the  speed 
and  incident  direction  ©(  ^ (|)c  striking  unit  surface  area 
per  unit  time.  Thus  Eq,  (4.6)  states  that  the  net  tangential 
momentum  transferred  from  the  surface  in  the  x direction  per  unit 
time  can  be  represented  by  the  integral  over  all  speeds  and  directions 
of  the  number  per  unit  time  of  molecules  incident  with  given  speed 
and  direction  multiplied  by  the  probability  that  these  molecules  will 
be  scattered  with  a given  speed  and  direction  multiplied  by  the  x 
component  of  that  speed  and  also  by  the  molecular  mass. 

Now  consider  how  (4*6)  may  be  specialized  by  application  of  the 
assumptions.  By  the  terms  of  assumption  (E-l)  Gr  must  be  determined 
by  the  same  speed  distribution  in  a particular  direction  as  that 
holding  for  molecules  issuing  from  a gas  in  equilibrium  at  the 
temperature  of  the  surface.  However,  this  implies  that  this  dis- 
tribution in  speed  is  independent  of  the  probability  that  an  issuing 
molecule  will  take  a given  direction  as  determined  from  the  measured 
flux  distributions  of  scattering  from  a real  surface.  The  issuing 
distribution  in  speed  is  also  assumed  by  (E-l)  to  be  independent  of 
the  incident  speed  distribution.  Under  the  terns  of  assumption  B 
the  probability  that  an  issuing  molecule  will  take  a given  direction 
is  also  assumed  to  be  independent  of  the  incident  speed.  Thus  the 
probability  distribution  function  P may  be  rewritten  as 

"Pi  Oc  (j)i 


It  will  be  shown  at  a later  point  in  this  section  that  P^  is  inde- 
pendent of  the  issuing  direction,  i0e.,  P0  = P0  (C  )<, 

3 j r 

It  will  be  recalled  that  the  representative  scattering  shells 
(see  section  3.3  and  Figure  42)  give  scattering  information  with 
respect  to  a primed  coordinate  system  with  the  X;  coordinate  axis 
oriented  along  the  incident  ray.  From  these  shells  a probability 
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distribution  function  of  the  form  P ( ©1,0'^')  can  be 
computed.  The  function  P ( ©i,©',  (j)')  is  assumed  (assumption  A) 
to  be  independent  of  the  orientation  of  X*  axis  in  the  scattering 
surface  and  hence  independent  of  (f)  l . A change  of  variables  can 
now  be  introduced.  Letting  (J)'  = and  O'  = ©r  ^ "P*. 

can  be  rewritten  as 

= ■PiCet,er,4>r+<t>i>  -?*  (©:.©', $) 


and  equation  (4.6)  becomes 

x sin*-  ©'  Sim  ©l  (4,7) 

In  which  the  angular  integrations  are  now  to  be  carried  out  over 
the  variables  4>*  • 

The  probability  distribution  Pg  may  be  expressed  in  terms  of, 
the  measured  flux  distributions  as 


(©;,©;<*>•)  = A(e,-.e‘,( !>') 


(4.8) 


A(©;  ,©',((>' )s\HG'dQ‘  d Cf>‘ 


where  Ate;©^)  is  the  beam  signal  amplitude  in  arbitrary  units, 
and  hence  is  proportional  to  the  molecule  number  flux  through  unit 
area  of  the  unit  scattering  hemisphere  in  the  direction  Q1 , (f)‘ 
i.e.,  At©^©'^'  ) = K N(  Oc.^©*,  (J)'  ) where  K is  a constant  of  pro- 
portionality. 


The  terms  in  the  primed  coordinates  may  be  collected  after  ex- 
pansion of  cos  (4>'“  (J>)  • The  integral  over  the  primed  coordinates  may 
be  written 


/ / ?2(&;#©)(j)'KtOS<t>'<-OJ(|>;  4 SlM^'uwfy)  SIM2©1 

Jo 
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(4.9) 
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It  should  be  noted  that  the  scattering  shells  are  symmetric  about  the 
X/  axis,  i.e.,  A ) = 4(9^0,-$),  an  even  function  in  (j)( 
about  * 0.  Hence  the  portion  of  the  integral  (4.9) 

/ I? 

J J SIN  SIN  (j>£ 


must  have  the  value  zero. 


The  remainder  of  (4.9)  may  be  written  as 


The  integrals  of  Cjy  , (defined  by  equation  (4.10)  must  be 
evaluated  numerically.  Details  of  the  numerical  integration  are  given 
in  Appendix  B and  a plot  of  the  results  is  shown  in  Figure  44. 


The  evaluation  of  (Cy)  must  be  considered  more  fully  at  this 
point.  Consider  a coordinate  system  \/Yj,  5 immersed  in  a body  of 
gas  at  rest  having  molecule  number  density  nf.  The  number  of  mole- 
cules per  unit  time  crossing  unit  area  of  the  plane  into  solid  angle 
(jlY.  making  angle  Ip  with  the  £ axis  and  having  speed  lying  within 
the  range  Cf  to  + d Cf  is  given  by 


Then  the  probability  per  unit  solid  angle  that  a molecule  will  have 
a speed  Cf  lying  within  unit  speed  range  will  be 

-Rlcv 

PU^  « Aa  e ef  cosi|) 

/ “+  A+  cosy 

If  this  fictitious  gas  is  assumed  to  be  thermal  equilibrium  with 
the  surface  then 

?U(.i=P3(cr)  = e'  r r c? 


'-a.vv 
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An  integration  of  P(Cr)  Cr  d over  all  speeds  may  be  performed  im- 
mediately such  that 


/ ?(Cr)  iCv  Cr-  37fT 

Jo  4 Pr 


(4.11) 


Equation  (4.7)  may  now  be  rewritten  as 

6r  » Am  Nl(Ci,0i/M  Cr  grC0Jcoi<t )i  SIN 
-4  Ci  (0i  , 

x4Cid0;  <4 (J) t 


(4.12) 


Equation  (4.12)  may  be  looked  at  in  a very  direct  and  intuitive 
fashion.  A review  of  the  integral  of  Eq.  (4.10)  defining  g ( 0;  ) 
will  show  that  the  quantity  wvCrQ»-(.0c)  might  be  thought  of  as  being 
the  magnitude  of  the  mean  X*  component  of  tangential  momentum  per 
molecule  transferred  from  unit  area  of  the  surface  per  unit  time  in 
the  primed  coordinate  system  by  molecules  incident  at  angle  Qi  • 
Multiplication  by  COS  ^ » » the  angle  made  by  the  X' 

axis  with  the  x axis  of  the  original  coordinate  system,  gives  the 
value  of  the  aforementioned  momentum  in  the  x direction  of  the  ori- 
ginal coordinate  system.  Multiplication  of  this  quantity  by  the  number 
of  molecules  incident  per  unit  time  of  all  speeds  from  all  directions 
and  an  integration  overall  vairables  yields  the  total  magnitude  of 
tangential  momentum  transferred  from  the  surface. 

To  return  to  the  original  argument,  if  the  integrand  of  equation 
(4.2)  is  put  into  polar  coordinate  form,  there  results,  as  the 
number  of  molecules  incident  on  the  surface  per  unit  time  from  di- 
rection 6^  and  (j)j  in  the  intervals  and  having  the  speed 
in  the  range  d C^,  the  expression 

- rv\  A C?5IN0*,COS©£ 

X e del  d©i, 


which  reduces  to 
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NUCi,0i,<t>i)  d(j>i  =■  /v\  A cf  Sin ©<.  c os©c 

_ [cl  - <2Ci.l)  siNJ0i.  coi(j);  + Ul  ^ 

* £ dcL  «i  ©i.  a<|>i. 


Then 


01  cos  <^; 


(4.U) 


Performing  the  indicated  integrations  over  and  (j);  we  have, 

l*  $*• 

Gy  ( 6;. ) =•  3YTT  0 * cos  ©i  & i ki  &i  a-  (©() 
lC>0r0  3 

* [Uk+2k3)I.^)+Um3)I,(£)] 

(4.15) 


where  b = s sin0{  and  IQ  and  1^  are  modified  Bessel  functions 
and  s(s  = |r)  i«  the  molecular  speed  ratio. 


The  details  of  the  integration  above  are  given  in  Appendix 

A and  — — ]/^  Is  plotted  as  a function  of  Qi 

for  s = 0.1,  1.0,  2,0,  and  10.0  as  given  in  Figure  43. 


The  final  integration  over  0{  must  be  performed  numerically  of 
course.  One  can  observe  immediately,  however,  that  for  gp  ( Q{  ) 

- 0 (diffuse  reflection)  Gr  will  equal  0 and  f * 1. 

Thus,  subject  to  the  limitations  of  the  data  and  of  the  assump- 
tions, f can  be  calculated  for  any  molecular  speed  ratio,  s,  from 
the  direct  measurement  of  the  scattered  flux.  Such  calculations  were 
carried  out  for  the  case  of  the  unpolished  glass  surface  for  values 
of  s = 0.1,  1.0,  2.0,  and  10.0  and  /3r  " 45.  A 

plot  of  the  resultant  findings  showing  the  behavior  of  f as  a function 
of  s is  shown  in  Figure  45. 


4.3  Calculation  of  f under  Terms  of  Assumption  E-2 
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The  main  work  of  this  section  is  again  the  calculation  of 
the  magnitude  of  the  tangential  momentum  transferred  from  the  gas 
to  the  surface  under  the  terms  of  assumption  E-2,  With  the  exception 
of  this  latter  assumption,  all  assumption^  are  as  outlined  in  section 
4*1*  The  same  flow  situation  is  also  assumed.  will  remain  as 
calculated  in  Eq.  4.5* 

Consider  now  again  how  the  probability  distribution  function  of 
Eq.  (4*6)  might  be  appropriately  specialized.  According  to  the  terms 
of  assumption  E-2  two  alternate  classes  of  gas-surface  interaction 
occur.  The  scattered  molecules  may  be  either  diffusely  scattered 
according  to  the  cosine  distribution  with  speed  distribution  as  in 
assumption  E-l  or  they  may  be  specularly  scattered  from  tilted  surface 
elements.  The  total  probability  distribution  may  be  represented  as 
the  sum  of  the  probability  distributions  for  these  two  alternate 
events.  Eq.  4.6  may  be  rewritten  as  follows* 


— ) }P|DUtpi^i,Cr,©r^Cr 

6{. 


+ 

x 


©r  (^r)GrJcOS  (j)r  SIN’-Gr 
d C‘i  d © t <1  (,  <jCf~  i 0|r  <1 


Sin  0t 
(4.16) 


The  subscript  D identifies  the  probability  function  associated  with 
the  diffusely  scattered  molecules  while  the  subscript  S identifies 
that  associated  with  the  specularly  scattered  molecules. 

By  the  same  arguments  as  were  employed  in  section  4.2, 
T>t0(Gi,es<K,cr,©r,Vl  can  be  specialized  to  the  form  T>lDce;e;<t>) 
kTjoCC*.)  .here  C|) ' = (J)r  + (|)L  and  0 ' - 0r  as  before. 

The  probability  distribution  function  P(i  (Q0  t<^(r>©r(|>r)  can  be  specialized 
immediately  to  the  form  Pzs(.9i>0$'  )■  since  the  angular  flux  dis- 
tribution is  assumed  to  be  independent  of  (assumption  B)  and  for 
these  molecules  Cr  can  have  only  the  value  C^.  Thus  equation  4.16 
can  be  rewritten  as 

Gr^  — Ni.CC i C©i.©;4),)1:>3p(Cr)Ct. 

+ ?*s(*«,©W'>Ci]cosCf-fe)  SIN1©’ 

x iiw  Qi  dc0  d (();  dCr  ©*  d 


(4.17) 
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Recalling  again  assumption  (E-2),  it  can  be  seen  that  A (j)'  ), 
the  beam  signal  amplitude  ( A (9;,©^  4)')  = K N(0t,O',(j)‘)  as  before), 
can  be  written 

A(e„,e',4>‘)  = adcos0'  + As  i <j>') 

where  Ap  is  a constant  representing  the  diameter  of  the  assumed 
diffuse  scattering  sphere  and  Ag  (0^0^  (j)')  are  the  amplitude  dif- 
ferences between  this  sphere  and  the  total  amplitudes  of  the  specular 
bulges.  From  the  data  then,  can  be  found  as  indicated 

below: 


-PzdO^')  Ad  cos  G' 

A (9^0*,  4)')  sim 

while  P2g(  0^0,^* ) becomes 


(4.18) 


V25(Q^$')= A&jgkJ^l (4.19) 

J J A (0t©',(j)'j  $jkj  Q'Jo'  d<(> 

o s <J 

Extracting  the  integrals  over  the  primed  coordinates  which  apply  to 
P2g  from  4«17,  we  have 


which  has  the  value  zero. 


The  remainder  of  the  integration  over  the  primed  coordinates 
may  be  written 


(4.21) 


51 


For  the  same  reasons  of  symmetry  as  before,  the  integration  over  the 
second  section  of  the  numerator  (l.e.  over  the  section 

S IN  <!>'***  0.)  makes  no  contribution  to  the  total 
values  of  the  integral.  Further,  we  can  add  A^  cos  0'  to 
AsCei,e',4>'  ) without  changing  the  value  of  the  integral.  Thus  the 
expression  4.21  can  be  re-expressed  as 

sXTt 

cos  ily  / A(&i  Q'CftlcosC^su^e^Q'cl  0 

(4>22) 

= COS  4>t  Cjrl©0 

The  <Jy/0{)  as  defined  here  is  the  same  as  that  defined  ty  Eq.  4.10 
of  the  previous  section. 

Equation  4-17  may  now  be  rewritten  as 


If  now  the  function  is  expressed  in  polar  coordi- 

nate form  as  in  equation  4.13,  and  if  the  integrand  of  4.23  is  multi- 
plied ty  — -■!  — §v — G r.  becomes 

py*'*  Act  © 

^ x.cos  Si^0;cos0;J  dCi,a©;  ^0* 


(4.24) 


or 


Gri=/  gir^Qj)  (4.25) 

/ft  Si  N 01 

where  ( Qi  ) is  the  integral  over  Ci  and  <$>;,  of  the  quan- 
tity within  the  brackets  of  Eq.  (4.24).  Although  this  grouping  might 
appear  rather  arbitrary  it  is  justified  by  the  fact  that  G|w(0^  has 
an  interesting  physical  interpretation  since  the  function  represents 
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the  net  x component  of  momentum  incident  on  the  surface  per  unit 
time,  unit  area  and  unit  polar  angle  at  angle  di  resulting  from 
molecules  of  all  speeds  coming  from  all  azimuth  angles* 

The  result  of  integration  over  and  (J)|  can  be  written  in 
the  form 


G; 


/YWrt'xVwv  cos  Ql  S 

X [(  J’"+AU3  + 2b)] 


(4*26) 


where 


=»  S b fe  * b — s jin®;  . s *»  ii_ 

aw 


and 


O A 1 , . 

O — / COS  (p*L  ERF  ( bcos  Q)t  ) e <J<fi 

^ (4.27) 


Details  of  the  integration  may  be  found  in  Appendix  C and  a 
plot  of  S“  may  be  found  in  Figure  46  giving  J*'  for  values  of  b from 
0 to  10.  Also  1183  teen  Plotted  as  a 

'Vvn  /v\  VwvT 

function  of  0£  for  s = 0.1,  1.0,  2.0,  and  10.0;  see  Figure  47. 

Thus  again,  subject  to  the  limitations  of  the  assumptions,  f 

can  be  calculated  for  any  molecular  speed  ratio.  Such  calculations 

were  performed  for  the  unpolished  glass  surface  for  s = 0.1,  1.0, 

2.0/  and  10.0  and  a plot  of  the  results  is  shown  in  Figure  45«  The 

final  integration  over  0 1 in  the  evaluation  of  G -had  to  be  per- 

r<c 

formed  numerically,  of  course. 

4.4  Calculation  of  f for  Pure  Specular  Reflection 

ai..Grazj^.Jprtflen£g 

A consideration  of  Figure  47  will  show  that  molecules  incident 
at  angles  close  to  grazing,  play  an  increasingly  important  part  in 
the  momentum  transfer  to  the  surface  as  the  molecular  speed  ratio 
increases.  It  has  already  been  argued  that  at  some  incident  angle 
very  close  to  grazing  a condition  of  pure  specular  reflection  may 
set  in.  (i.e.,  the  reflected  molecule  has  direction  <{>£  andfT- 
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and  has  speed  Cr  * C^).  All  molecules  at  some  angle  S©  on  to  90° 
will  be  assumed  to  be  specularly  reflected  for  the  purpose  of  this 
speculative  calculation  while  all  molecules  incident  form  Ct  to  e. 
will  be  assumed  to  be  diffusely  reflected  according  to  the  cosine  law 
and  thus  contribute  nothing  to  the  momentum  transfer  from  the  surface, 
f for  this  situation  can  be  very  simply  formulated  as 

Gi.  (4.28) 

where  the  tangential  momentum  transferred  from  the  surface  is  equal 
in  magnitude  to  G»U0i)  integrated  from  ©4  to  . Calculations 

of  f were  performed  for  assuming  various  values  from  85°  to 

90#  and  for  s,  the  molecular  speed  ratio  taking  on  the  values  s = 0.1, 
1.0,  2.0  and  10.0.  The  results  of  these  calculations  are  plotted  in 
Figure  48. 

4.5  Discussion 

Figure  45  shows  the  behavior  of  f for  variations  of  the  para- 
meter s for  both  assumed  cases  of  molecular  interaction  E-l  and  E-2. 

The  diverging  behavior  of  these  plots  lends  plausibility  to  the 
intuition  which  had  originally  led  to  assumptions  E-l  and  E-2;  namely 
that  these  assumptions  were  bracketing  assumptions  and  that  more  com- 
plete knowledge  of  the  scattering  probability  distribution  function 
would  probably  lead  to  values  of  f for  this  surface  lying  in  be- 
tween these  extremes. 

It  is  possible  that  the  relative  validity  of  assumptionsE-1  and 
E-2  leading  to  the  divergent  results  above  might  be  assessed  by  the 
performance  of  a rotating  cylinder  experiment  at  molecular  speed 
ratios  of  2 or  more.  Unfortunately,  such  an  experiment  would  be 
difficult  to  perform  because  of  the  high  peripheral  speeds  Involved 
and  the  delicate  drag  measurements  necessary  to  resolve  the  rather 
small  differences  in  f.  A much  more  reasonable  approach  would  be, 
as  suggested  in  section  3.4»  to  introduce  a velocity  selector  into 
the  path  of  the  scattered  molecules  in  a molecular  beam  experiment 
and  measure  directly  the  reflected  speed  distributions. 
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It  is  interesting  to  note  that  neither  calculation  based  on  the 
alternate  assumptions  (E-l)  and  (E-2)  yielded  values  of  f nearly  loir 
enough  for  agreement  with  the  reported  value  of  f * 0.89  for  air  on 
glass.  Such  a value  is  clearly  beyond  the  limit  of  error  of  these 
present  experiments  since  it  could  only  result  from  much  larger  for- 
ward distortions  of  the  scattering  shells.  Millikan,  in  arriving  at 
his  value  of  f (above)  re-examined  the  earlier  work  of  Knudsen  and 
of  Warburg  on  slip  in  glass  capillaries.  Warburg  worked  with  air  on 
glass  while  Knudsen  confined  his  investigations  to  H2,  02  and  C02. 

It  is  barely  possible  that  the  somewhat  higher  pressures  at  which  these 
early  experiments  were  performed  has  some  effect  on  the  adsorbed  gas 
layer  thickness,  which  in  turn  may  have  allowed  more  "specular*  scat- 
tering. It  is  also  likely  that  some  errors  could  be  expected  as  a re- 
sult of  using  the  findings  of  experiments  on  gases  other  than  air  to 
make  predictions  for  the  air-glass  surface  Interaction.  However,  no 
explanation  has  been  found  which  properly  accounts  for  the  discrepancy 
between  Millikan's  results  and  that  of  this  report. 

The  results  of  the  calculation  of  section  4.4  are  also  quite 
interesting  since  they  serve  to  show  the  importance  of  the  gas-surface 
behavior  at  near  glancing  angles  of  incidence  in  determining  the  over- 
all momentum  transfer.  These  results  also  serve  to  show  the  impor- 
tance of  extending  the  molecular  beam  scattering  experiments  to  in- 
clude measurements  at  angles  extremely  close  to  glancing.  It  need 
hardly  be  stated  that  an  engineering  surface  which  exhibited  some 
considerable  degree  of  specular  reflection  at  a few  degress  from 
glancing  incidence  would  find  immediate  application  in  high  altitude, 
high  speed  aerodynamics. 

It  is  hoped  that  the  calculations  of  this  section  will  suggest 
a new  approach  to  the  problems  of  evaluating  the  boundary  conditions 
for  rarefied  gas  flows.  It  seems  possible  that  the  old  concepts  of 
a constant  f and  a constant  oc  , the  energy  transfer  counterpart  of 
f,  might  profitably  be  supplanted  by  the  more  general  scattering  pro- 
bability function. 

5.0  SUMMARY 

A molecule  r beam  apparatus  was  used  to  investigate  the  scatter- 
ing of  air  and  N2  molecules  from  surfaces  of  polished  low  carbon  steel, 
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etched  low  carbon  steely  polished  aluminum,  unpolished  window  glass,  and 
polished  window  glass.  A beam  of  molecules  was  directed  against  the  sur- 
faces at  representative  angles  of  incidence.  Polar  flux  distributions 
of  the  scattered  molecules  were  measured  by  a movable  detector.  An  ioni- 
zation gage  beam  detector  system  was  devised  to  accomplish  this  measure- 
ment. 

Spatial  polar  plots  of  these  measured  flux  distributions  were 
found  to  approximate  closely  the  form  of  the  cosine  scattering  distribu- 
tion for  the  cases  of  the  steel  and  aluminum  surfaces.  Similar  plots  in 
the  case  of  the  glass  surfaces  were  found  to  have  well  defined  but  small 
bulging  deviations from  the  cosine  scattering  shape  on  the  side  away  from 
the  incident  beam.  These  findings  did  not  serve  to  distinguish  the  degree 
of  elasticity  of  the  interaction  since  both  elastic  and  inelastic  inter- 
action models  would  serve  to  produce  pure  cosine  scattering,  and  only  a 
small  degree  of  elasticity  would  suffice  in  some  models  to  account  for 
the  observed  deviation  from  cosine  scattering.  Among  others,  a modified 
Maxwell  model  would  seem  to  satisfy  the  available  data.  In  this,  a 
fraction  of  the  incident  molecules  undergo  a diffuse  scattering,  being 
more  or  less  accomodated  in  energy  to  the  surface,  while  the  remainder  are 
specularly  reflected  from  tilted  surface  elements. 

In  order  to  illustrate  one  aerodynamic  application  of  scattering 
data  obtained  by  molecular  beam  techniques,  a procedure  was  devised  for 
the  computation  of  the  coefficient  of  specular  reflection,  f,  and  values 
of  this  coefficient  were  calculated  for  the  case  of  an  unpolished  glass 
surface  oriented  tangentially  to  the  direction  of  mass  motion  in  a free 
molecule  flow.  These  values  were  calculated  as  a function  of  the  para- 
meter s,  the  molecular  speed  ratio.  If  all  scattered  molecules  are  as- 
sumed to  possess  the  distribution  in  speed  in  a particular  direction 
which  they  would  have  had  had  they  issued  from  a gas  in  equilibrium  at 
the  temperature  of  the  surface,  f takes  on  values  of  rj  0.97  as  s approaches 
0 to  /v  0.99  at  s = 10.  Alternatively,  if  the  modified  Maxwell  surface 
interaction  model  is  assumed,  f takes  on  values  of  rJ  0.97  as  s approaches 
0 to/\/  0.93  at  s 3 10.  These  values  do  not  agree  with  the  published 
value  of  f of  0.89  for  air  on  glass.  This  lack  of  agreement  could  not  be 
accounted  for  by  a consideration  of  experimental  errors. 

A speculative  calculation  was  performed  showing  the  effect  on  f 
for  molecular  speed  ratios  of  0.1,  1.0,  2.0,  and  10.0,  in  the  case  where 
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a condition  of  pure  specular  reflection  occurs  from  source  angle,  Q©  9 
close  to  grazing  incidence  on  to  90°  incidence,  f to  a shown  to  have  the 
value  *0.5  for  s ■ 10  and  ©o  = 85*  • 

Three  points  of  departure  for  further  experimental  work  were 
discussed,  A.  The  substitution  of  mercury  diffusion  pumps  for  the  pre- 
sent oil  diffusion  pumps  would  remove  one  important  source  of  surface  con- 
tamination. B.  Higher  resolution  experiments  using  a slit-form  beam 
geometry  would  yield  interesting  data  on  the  behavior  of  the  scattering 
at  incident  angles  close  to  grazing.  C.  A velocity  selector  introduced 
into  the  scattered  beam  would  afford  a means  of  determining  the  relative 
validity  of  the  assumptions  leading  to  the  computed  values  of  f.  This 
device  would  also  afford  a means  of  studying  in  detail  the  energy  exchange 
between  gas  and  surface.  This  latter  work  would  provide  a valuable  basis 
for  a better  physical  tinderstanding  of  the  basic  scattering  mechanisms. 
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ieebhmla 

To  integrate  over  Ci  and 


(A-l) 


i *1  i 


let  "P -<2.Us\N©CO,$<^ 

and  ££.=  LLmn£coc(|>  * bCO 

2-  Vim 


Making  the  above  substitution,  completing  the  square  in  the  exponent, 
and  letting  /yvn/v\_  3/tFA  cos  O SIM  Q frt- (£)  =s  N4 

4 (Sr 


there  results 


P c-  1,1  Jc 


Integration  over  C yields 

%Tf 


° - - --  ^ pi 

Gr(.S)(i6  = M<ie/c.t>i(()e'B  ('U  'tI  LgE  |\  +e.RF(&&)j 

+ + w(¥)l}]4  <«> 


Multiplying  out  and  collecting  terms  it  was  found  that  the  integrals  over 
(|>  not  equal  to  zero  ares 

_ /lTT  i . 

3Vff  /coslq)e 

1 P U ,x  - 

= b [l„t£)  + !/§)] 


t P+ 


U-3) 


and 
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^rr  ^tos’-cj) 
cos*(y  6 49 


(A-4) 
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I0(f  jandl^-j ) are  the  zeroth  and  first  order  modified  Bessel  functions 
of  the  first  kind. 

A final  collection  of  terms  yields  the  result  of  Eq.  4.15: 

joL_  i2- 

Gr(^t3  s 3YfT  6 1 cos  ©l  simQl 

I6(*rf3  9 

X Ijj  l>+2t3)Xo^)  + Uo+2b3)I,^  )j 
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APPENDIX  B 

In  order  to  perform  the  indicated  integration  of  Eq.  4.11  l.e. 

¥ 

A ) cos  § <J©  d(|) 

and 


the  indicated  integrations  must  be  replaced  by  summations  and  the  sum- 
mation carried  out  over  the  surface  of  a unit  hemisphere.  Thus  the  flux 
through  a surface  element  is  approximated  by  A^lOi^d)' ) A Q.^  where 

Aa^  a = a Tr  (yN>r4;x 

<v\-0 

A simple  integration  scheme  was  devised  in  which  the  total 
unit  hemisphere  was  divided  into  414  elementary  surfaces  of  approximately 
equal  area.  Since  the  scattering  distributions  were  assumed  symmetrical 
about  the  incident  plane,  just  half  of  this  number  were  employed  in  the 
actual  integration.  Since  all  scattering  shells  were  regarded  as  spherical 
on  the  incident  ray  side,  this  section  of  the  integration  was  performed 
only  once.  For  a representation  of  the  division  of  the  unit  hemisphere 
surface  see  Figure  49. 

The  were  taken  from  the  representations  of  the  scattering 

shells  at  the  centerpoint  of  each  elementary  area.  The  resultant  ^ are 
illustrated  in  Figure  44. 

A table  of  the  AG.* of  this  integration  scheme  follows: 

BOUNDING  POLAR  ANGIES  NUMBER  OF  ELEMENTARY  AREA  OF 

OF  SPHERICAL  ZONE  SURFACES  IN  SPHERICAL  ZONE  SURFACE  ELEMENT 

90° 

85° 

80° 

75° 

70° 

65° 

60° 

55° 


85° 

36 

0.015212 

80° 

36 

0.015048 

75° 

36 

0.014866 

70° 

36 

0.014520 

65° 

36 

0.014068 

o 

O 

vO 

36 

0.013503 

55“ 

36 

0.012843 

50® 

36 

0.012083 
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BOUNDING  POLAR  ANGLES 
OP  SPHERICAL  ZONE 

NUMBER  OF  ELEMENTARY 
SURFACES  IN  SPHERICAL  ZONE 

AREA  OF 

SURFACE  ELEMENT 

50?  - 40? 

36 

0.024586 

40?  - 30° 

36 

0.017455 

o 

O 

OJ 

1 

o 

o 

30 

0.015425 

20°  - 10° 

18 

0.015752 

10°  - 0° 

6 

0.015909 

65 


If^T 

v e" 


APPENDIX  C 

To  evaluate  the  net  Incident  tangential  momentum  Gi<w)  as 
a function  of  0;  where 

Gi  (0,\  =t  nnr><v\  Q3  S>m'uQlCO$Q,‘.  / / C*CO$()>i 

ir^i  A Jo 

£c*-2C.i,UstN0uCoi  Cj>i  + 

(Prom  Eqs.  4.24  and  4.25) 

let  T * Umn  0cos  (|j 
and  P(3  = U sumQcqsC})  «. 

2 V!w\ 
as  in  Appendix  A. 

Making  the  above  substitution,  completing  the  square  in  the 
exponent  and  letting  # 

/W\  /W  S|Wl©COi  © ss.  M 


ir3A 


there  results 


Gl(e, . M fT^ .'^V e ^ l)Vd<|> 


Integration  over  C yields 

, -e’Ur-f) 


G;(e>=  cos<j>e 

(e  [i  + «,V 


1 

+• 


1 (Q-2) 


Multiplying  out  and  collecting  terms  the  integrals  over  <})  not  equal 
to  zero  were  found  to  bes 


k3  e~  / <j)d(b  g 5 tt  b3  6 
2 a*  Jo  SO* 


-zir 


- 5 


O 

.2  r* n 


sb _&7  cosld)d(b  = r Tr  b e~ 

Gy  Jo 

-SfiLeT4’*1 y/cos  C(>  ERF(bcosd))e-^  ^dCj) 

8Gff  ° J* 


=.  2W  e 

0 e5- 


(C-3) 

(C-4) 


(C-5) 
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ayf5’^/  cos3<J>s»F(Ww<|>le-k'il"’l*,4$ 

= m'e'1’*11  p_f(2bM)lJ0’ 

Isf  U*>  L J J (C-6) 

nr  b+  e"  ^7o  4>  EKF(bcos<(>)e'b'SIN^d<t> 

-^“VE1  +we'l‘[4-AJ)  (0-7) 

Integrations  (C-6)  and  (C-7)  proceed  by  parts  from  (C-5),  A 
collection  and  cancellation  of  terms  yields 

GlCe)®J^e“Slpfreb\y^4  b'j  + fr  b3  +2TrbJ 

If  now  <r  is  defined  such  that  S’-  h^S 

2ft T 

a final  substitution  yields  the  results  of  Eq.  4*26: 

Gl(Q;)=  nrv\/y\  Vv^  COS  Ol  0 * [ (^"^  + 0(  k3  + 2.  t> jH  (C-8) 

Ytt  ** 

There  remains  the  task  of  evaluating  ■S  or  This  integral  did  not 

appear  reducible  to  previously  tabulated  functions.  A numerical  inte- 
gration was  accomplished  as  shown  below: 

COS  (J)  EK F (t>COi<jp)e~^  ^d<|) 

let  X = b CO 3 (t)  «,  d(b  * — <j  * 

^ Yb*-XT 


The  appropriate  transformation  of  the  limits  and  an  inte- 
gration by  parts  makes  possible  the  expression  of  S as 

S = 2e'wWf  b-1  -yL  / Vbv-Xl  En f (x ) e*' d* 
b L 

applying  the  definition  of  >0  there  results 


cT  - bl+_4_  ERF(X)e*  AA 

Wl 


(C-9) 


In  this  form  S'  was  computed  for  values  of  b sufficient  for  the 
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further  computation  of 

r-ft  ryf  'VW*' 

of  O and  logo  follows  and 

i 


A table  of  the  computed  values 


/VW(V\  Vy 

a plot  of  these  values  is  given  in  Fig.  46. 


<S" 


LOfa 


S' 


0 

0.0000 

- 00 

0.005 

0.00249 

-2 . 6038 

0.1 

0.01005 

-1.9978 

0.15 

0.02276 

-1.6428 

0.2 

0.04081 

-1.3892 

0.25 

0.06449 

-1.1905 

0.3 

0.09417 

-1.0261 

0.35 

0,13032 

-0.8850 

0.4 

0.17351 

-0.7607 

0.45 

0.22446 

-0.6489 

0.5 

0.28402 

-0.5466 

0.55 

0.35322 

-0.4520 

0.6 

0.43330 

-0.3632 

0.65 

0.52571 

-0.2792 

0.7 

0.63219 

-0.1992 

0.75 

0.75480 

-0.1222 

0.8 

0.89598 

-0.0477 

0.85 

1.05864 

0.0248 

0.9 

1.24023 

0.0956 

0.95 

1.46284 

0.1652 

1.0 

1.71331 

0.2338 

1.25 

3.772 

0.5766 

1.5 

8.478 

0.9283 

1.75 

20.37 

1.309 

2.0 

53.58 

1.729 

4.0 

8.891  x 106 

6.949 

10.0 

2.69  x 1043 

43.43 

TABLE  I 

COEFFICIENT  OF  SPECULAR  REFLECTION  f 
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GAS 

SURFACE 

f 

REFERENCE 

air 

machined  brass 

1.00 

11 

co2 

machined  brass 

1.00 

11 

air 

old  shellac 

1.00 

11 

co2 

old  shellac 

1.00 

11 

air 

Hg 

1.00 

11  ■ 

air 

oil 

0.90 

11 

CM 

o 

o 

oil 

0.92 

11 

H2 

oil 

0.93 

11 

air 

glass 

0.89 

11 

He 

oil 

0.87 

11 

air 

fresh  shellac 

0.79 

11 

air 

ig2° 

0.98 

12 

He 

ig20 

1.00 

12 

H2 

4g20 

1.00 

12 

°2 

4g2° 

0.99 

12 

air 

oil  on  machined 
aluminum 

0.90 

13 

TABIE  II 

THERMAL  ACC CHMQDATI ON  COEFFICIENTS  «C 
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GAS 

SURFACE 

CK 

REFERENCE 

H2 

bri$it  Pt. 

0.32 

14 

H2 

black  Pt. 

0.74 

14 

°2 

bright  Pt. 

0.81 

14 

°2 

black  Pt. 

0.93 

14 

N2 

Pt 

0.50 

15 

n2 

tungsten 

0.35 

15 

air 

flat  laquer  on  bronze 

0.88  - 0.89 

16 

air 

polished  bronze 

0.91  - 0.94 

16 

air 

machined  bronze 

0.89  - 0.93 

16 

air 

etched  bronze 

0.93  - 0.95 

16 

air 

polished  cast  iron 

0.87  - 0.93 

16 

air 

machined  cast  iron 

0.87  - 0.88 

16 

air 

etched  cast  iron 

0.89  - 0.96 

16 

air 

polished  aluminum 

0.87  - 0.95 

16 

air 

machined  aluminum 

0.95  - 0.97 

16 

air 

etched  aluminum 

0.89  - 0.97 

16 

He 

tungsten 

0.025  - 0.057 

17 

He 

nickel  not  flashed 

0.20 

17 

He 

nickel  flashed 

0.085 

17 

He 

tungsten  flashed 

0.17 

18 

0.12 

15 

He 

tungsten  not  flashed 

0.53 

18 

argon 

tungsten  flashed 

0.82 

18 

0.46 

18 

argon 

tungsten  not  flashed 

1.00 

18 

TABLE  III 

STM! ART  OF  MOLECULE-SURFACE  SCATTERING  INVESTIGATIONS 
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FIG.  I GEOMETRY  OF  SCATTERING  EXPERIMENT 
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SOURCE  REGION 
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RECORDER 


FIG.  3 MOLECULAR  BEAM  APPARATUS 
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FIG. 4 OPTICAL  BENCH  WITH  BEAM 
FORMING  AND  DETECTING 
EQUIPMENT 
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FIG. 5 SCHEMATIC  OF  SOURCE  GAS  AND  CALIBRATION  SYSTEM 
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FIG.6  SCHEMATIC  OF  BASIC  ROTATION  ELEMENTS 
OF  DETECTOR  MECHANISM 
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FIG.  8 IONIZATION  GAGE-MOLECULAR 
BEAM  DETECTOR 
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FIG.  10  BLOCK  DIAGRAM  OF  EMISSION  REGULATOR 
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FIG.  II  MOLECULAR  BEAM  ION  GAGE  REGULATOR 
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CHASSIS  CONNECTIONS: 
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FIG.  12  b 
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SPECIMEN  HOLDER 
METAL  SURFACES 


ETCHED  CRS 
SURFACE 
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SURFACE 


SPECIMEN  HOLDERS 
GLASS  SURFACES 


POLISHED  ALUM- 
INUM SURFACE 


UNPOLISHED  GLASS 
SURFACE 


POLISHED  GLASS 
SURFACE 


FIG.  16  TEST  SURFACES  & HOLDERS 
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FIG. 17  RAW  DATA  RECORD 

(FROM  RUN  OF  JULY  22,  1953  ) 
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FIG.  18  SCATTERING  PLOT-  POLISHED  STEEL  SURFACE 


FIG.I9  SCATTERING  PLOT- POLISHED  STEEL  SURFACE 
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FIG. 20  SCATTERING  PLOT-POLISHED  STEEL  SURFACE 
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INC  ANGLE  0.00° 

AZ  ANGLE  0.00° 
UNHEATED  SPECIMEN 
o-  MARCH  23, 1953 
HEATED  SPECIMEN 
□ -APRIL  7,  1953 
UNHEATED  SPECIMEN 
a - APRIL  8,  1953 


FIG.  22  SCATTERING  PLOT-ALUMINUM  SURFACE 

INC  ANGLE  25.11° 


FIG.  23  SCATTERING  PLOT- ALUMINUM  SURFACE 
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o-MARCH  23,  1953 
HEATED  SPECIMEN 
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FIG.  24  SCATTERING  PLOT-  POLISHED 
ALUMINUM  SURFACE 


FIG.25  SCATTERING  PLOT-  POLISHED 
ALUMINUM  SURFACE 
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INC  ANGLE  0.00° 


AZ  ANGLE  0.00° 
a MAY  16,  1953 

INC  ANGLE  0.00° 
AZ  ANGLE  19.57° 
° MAY  23,  1953 


FIG.  26  SCATTERING  PLOT  UNPOLISHED 
GLASS  SURFACE 
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° MAY  16,  1953 
o MAY  5,  1953 

INC  ANGLE  25.11° 
AZ  ANGLE  19.57° 
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FIG.  27  SCATTERING  PLOT  UN  POLISHED 
GLASS  SURFACE 
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INC  ANGLE  40.67° 
AZ  ANGLE  35.83° 
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FIG.  28  SCATTERING  PLOT  UNPOLISHED 
GLASS  SURFACE 


FIG.  29  SCATTERING  PLOT  UNPOLISHED 
GLASS  SURFACE 
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FIG  30  SCATTERING  PLOT  UNPOLISHED 
GLASS  SURFAFE 


FIG.  31  SCATTERING  PLOT  UNPOLISHED 
GLASS  SURFACE 
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FIG.  32  SCATTERING  PLOT  UNPOLISHED 
GLASS  SURFACE 
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FIG. 33  SCATTERING  PLOT- POLISHED 
GLASS  SURFACE 


FIG.  34  SCATTERING  PLOT  - POLISHED 
GLASS  SURFACE 
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FIG. 35  SCATTERING  PLOT- POLISHED 
GLASS  SURFACE 


INC.  ANGLE  56.25° 
AZ.  ANGLE  0.00° 

o-JULY  23,1953 
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FIG.  36  SCATTERING  PLOT- POLISHED 
GLASS  SURFACE 


HYD  6556-150-0 


FIG. 37  SCATTERING  PLOT- POLISHED 
GLASS  SURFACE 


INC.  ANGLE  79.62° 
AZ.  ANGLE  0.00° 
°-JULY  22, 1953 


FIG.  38  SCATTERING  PLOT- POLISHED 
GLASS  SURFACE 
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FIG.  39  SCATTERING  PLOT  - POLISHED 
GLASS  SURFACE 


FIG. 40  SCATTERING  PLOT- POLISHED 
GLASS  SURFACE 
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FIG.4I  COORDINATE  SYSTEM  - INCIDENT 
MOMENTUM  CALCULATIONS 
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FIG.  42  COORDINATE  SYSTEM  - 

ISSUING  MOMENTUM  CALCLATIONS 
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FIG.  45  f (COEFFICIENT  OF  SPECULAR  REFLECTION)  AS  A 
FUNCTION  OF  s (MOLECULAR  SPEED  RATIO)  FOR 
NITROGEN  ON  UNPOLISHED  GLASS 


HYD  6567-150*0 


FIG.  49  DIAGRAM  OF  UNIT  HEMISPHERE  AREA  DIVISION 
FOR  ISSUING  MOMENTUM  CALCULATIONS 


DISTRIBUTION  LIST  FOR  0 NR  TECHNICAL  REPORTS 

Director,  Aeronautical  Res,  NACA,  1724  P St.  NW,  Wash.  25,  D.  C.  (l) 

Commanding  General,  Air  Res.  and  Development  Command,  PO  Box  1395,  Baltimore  (2) 
Maryland,  Attn:  OSR(RDRF) 

Dr.  Morton  Alperin,  Western  Regional  Office,  Hdqrs,  Air  Res.  & Development  (l) 
Command,  PO  Box  2035,  Pasadena  2,  Calif. 

Director,  Ames  Aeronautical  Lab.,  NACA,  Moffett  Field,  Calif.  (l) 

Applied  Physics  Lab.,  Johns  Hopkins  Univ.,  8621  Georgia  Ave,  Silver  Sprg,  Md.(l) 

Armed  Forces  Special  Weapons  Project,  The  Pentagon,  Wash.,  25,  D.  C.  (l) 

Commanding  General,  AF  Cambridge  Res.  Cen.,  230  Albany  St,  Attn:  Dr.  H.  D.  (l) 
Edwards',  Geophysics  Div.,  Cambridge  39,  Mass. 

Commanding  Officer,  Arnold  Engr.  Devel.  Center,  Tullahoma,  Tenn.,  Attn:  (l) 

ACGT,  ACR  & AEMCL 

Mr.  C.  K.  Bauer,  Hermes  Proj.,  General  Electric  Co,  Schenectady  5,  N.Y.  (l) 

Dr.  Seymour  Bogdonoff,  James  Forrestal  Res.  Cen,  Princeton  Univ,  Princeton,  (l) 
N.  J. 

Prof.  W.  S.  Bradfield,  Dept,  of  Aeronautical  Engr.  Univ.  of  Minnesota  (l) 

Minneapolis,  Minn. 

Chief,  Bureau  of  Aeronautics,  Dept,  of  the  Navy,  Wash.,  25,  D.C.  - Attn:  (l) 

Aero  and  Hydro  Branch 

Chief,  Bureau  of  Aeronautics,  Dept,  of  the  Navy,  Wash.  25,  D.C.  - Attn:  (l) 

Research  Division 

Bureau  of  Aeronautics  Representative,  USN,  Central  District,  Wright-Patter-  (l) 
son  AFB,  Dayton,  Ohio,  Attn:  WC0WS-3,  James  G.  Law. 

Chief,  Bureau  of  Ordnance,  Dept,  of  the  Navy,  Wash.  25,  D.C.  - Attn:  Re9a  (l) 

Dr.  Francis  H.  Clauser,  Johns  Hopkins  Univ.  Dept,  of  Aeronautical  Engr.,  (l) 
Baltimore  18,  Md. 

Commanding  Officer  and  Director,  David  Taylor  Model  Basin,  Wash.  7,  D.C.  - (l) 

Attn:  Technical  Library 

Chief  of  Naval  Res.,  Dept,  of  the  Navy,  Wash.  25,  D.C.  - Attn:  Code  438  (3) 

Chief  of  Staff,  Hdqrs  USAF,  Directorate  of  Research  and  Development,  Wash.  (l) 

D.C.  - Attn:  AFDRD-AC2 

Chief,  Documents  Service  Center,  Armed  Services  Technical  Information  (2) 

Agency,  Knott  Bldg.,  Dayton  2,  Ohio  (1  reproducible  copy) 

Prof.  Fred  S.  Eastman,  Univ.  of  Wash.,  Dept,  of  Aero.  Engr.,  Seattle  5>  Wash.(l) 

Dr.  I.  Estermann,  Code  419,  Office  of  Naval  Research,  Wash.,  D.C.  (l) 

Prof.  K.  0.  Friedrichs,  New  York  Univ,  Inst,  for  Math,  and  Mechs.,  45  (l) 

Astor  Place,  New  York,  N.  Y. 

Director,  Jet  Propulsion  Lab.,  Calif,  In3t.  of  Tech.,  Pasadena,  Calif.  (l) 

Dr,  J.  M.  Jones,  Engr.  Res.  Inst.,  Univ.  of  Mich.,  Ann  Arbor,  Mich.  (l) 

Prof.  A.  Kantrowitz,  Graduate  School  of  Aeronautical  Engr.,  Cornell  Univ.  (l) 
Ithaca,  New  York. 

Dr.  J.  Kaplan,  Dept,  of  Physics,  Univ.  of  Calif,  at  Los  Angeles,  L.A.,  Calif. (l) 

Prof.  J.  H.  Keenan,  Mass.  Inst,  of  Tech.,  Mechanical  Engr.  Dept,  Cambridge  (l) 
Mass. 


Nov.  1953  DISTRIBUTION  LIST  FOR  ONR  TECHNICAL  REPORTS  . Page  2 

Prof.  M.  Z.  Krsywoblocki,  Univ.  of  Illinois,  Urbana,  Illinois  (l) 

Dr.  A.  M.  Kuethe,  Dept,  of  Aero.  Engr.,  Univ.  of  Mich,  Ann  Arbor,  Mich.  (l) 

Dr.  A.  R.  Kuhlthau,  Univ.  of  Virginia,  Dept,  of  Physics,  Charlottesville,  Va.(l) 
VIA:  Naval  Inspector  of  Ordnance,  Applied  Physics  Lab.,  Johns  Hopkins 
Univ.,  Silver  Spring,  Md. 

Director,  Langley  Memorial  Aeronautical  Lab.,  NACA,  Langley  Field,  Virginia  (l) 

Prof.  Lester  Lees,  Dept,  of  Aero.  Engr.,  Calif.  Inst,  of  Tech.,  Pasadena  (l) 

Director,  Lewis  Flight  Propulsion  Laboratory,  NACA,  Cleveland,  Ohio  (l) 

Prof.  H.  Leny,  Applied  Math.  Group,  Math.  Dept.,  Univ.  of  Calif.,  Berkeley  (l) 

Dr.  James  E.  Lipp,  The  Rand  Corp.,  1500  Fourth  St.,  Santa  Monica,  Calif.  (l) 

Prof.  J.  R.  Markham,  Mass.  Inst,  of  Tech.,  Cambridge  39,  Mass.  (l) 

Dr.  C.  B.  Millikan,  Guggenheim  Aeronautical  Lab.,  Calif.  Inst,  of  Tech.  (l) 

Pasadena,  California 

Dr.  Henry  T.  Naganatsu,  Director,  Hypersonic  Wind  Tunnel,  Calif.  Inst,  of  (l) 
Tech.,  Pasadena  4,  Calif. 

Commanding  Officer,  U.S.  Naval  Ordnance  Test  Station,  Inyokern,  China  Lake  (l) 
Calif.  - Attn:  Technical  Libraiy 

Director,  Naval  Research  Lab.,  Wash.  20,  D,  C.  - Attn:  Code  2020  (9) 

Commander,  Naval  Ordnance  Lab.,  White  Oak,  Silver  Spring  19,  Md,  Attn:  (l) 

Aeroballistic  Res.  Dept. 

Director,  Office  of  Naval  Research,  346  Broadway,  New  York  13,  N.  I.  (l) 

Director,  Office  of  Naval  Research,  844  North  Rush  St.,  Chicago  11,  111.  (l) 

Director,  Office  of  Naval  Research,  1000  Geary  St.,  San  Francisco  9,  Calif.  (2) 

Director,  Office  of  Naval  Research,  1030  East  Green  St.,  Pasadena  1,  Calif,  (l) 

Officer  in  Charge,  Office  of  Naval  Research,  Navy  #100,  FPO,  New  York,  N.Y.  (3) 

Radiation  Lab.,  Univ.  of  Calif.,  Berkeley  4,  Calif.  - Attn:  R.  Wakerling  (l) 

Chairman,  Research  and  Development  Board,  The  Pentagon,  Wash.  25,  D.C.  (l) 

Attn:  Information  Requirement  Board,  Room  3D -107 5 

Mr.  R.  H.  Shick,  Consolidated  Vultee  Aircraft  Corp.,  San  Diego  12,  Calif.  (l) 

Mr.  Joseph  Sternberg,  Ballistics  Res.  Lab.,  Ordnance  Dept,  Aberdeen,  Md.  (l) 

Sverdr  j & Parcel,  Inc.,  Consulting  Engrs.,  Syndicate  Trust  Bldg,  St.  Louis  (l) 

Dr.  Theodore  Theodorsen,  Univ.  of  Maryland,  Bldg.  J,  Rm.  314,  College  Pk,  Md.(l) 

Dr.  F.  L.  Wattendorf,  A.G.A.R.D.  of  NATO,  APO  58,  c/o  PM,  New  York  (l) 

Western  Coordination  Office,  NACA,  7660  Beverly  Blvd.  L.  A.,  36,  Calif.  (l) 

Dr.  E.  P.  Williams,  Missiles  Div.,  Rand  Corp.,  1500  Fourth  St.,  Santa  (l) 

Monica,  California 

Wright  Air  Development  Center,  Wright-Patterson  Air  Force  Base,  Dayton,  Ohio  (2) 
Attn:  WCRRS-3  (l)  Attn:  MCIDXE  (l) 

Dr.  J.  H.  Marchant,  Div.  of  Engr.,  Brown  Univ.,  Providence  12,  Rhode  Island  (l) 


